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1 .0  INTRODUCTION 


1.1  GENERAL  CONCEPT  OF  COUNTERDOPED  DETECTORS 

The  development  of  a  counterdoped  infrared  detector  is  due  t( 
Elliott.^  ^  A  background  limited  infrared  photoconductor  (BLIP; 
which  operates  at  high  temperature  (80K)  is  made  possible  using  the 
technique  described  by  Elliott  et  al.  A  normal  infrared  detector  wi  1 ' 
exhibit  a  very  large  capture  cross-section  of  an  ionized  acceptor  for  i 
hole.  This  large  capture  cross-section  (attraction  of  carriers  anc 
photoionized  centers)  leads  to  short  carrier  lifetimes  which  result  ir 
low  operating  temperatures  for  the  detector.  The  counterdoped  syster 
shows  that  holes  are  captured  by  donors  which  are  neutral  or  positively 
charged;  the  capture  cross-section  in  this  case  is  very  small  because  th( 
photoionized  center  is  neutralized  in  the  photoabsorption.  Thi; 
resulting  lower  recombination  cross-section  increases  the  carrier 
lifetime  permitting  development  of  a  BLIP  detector  which  will  operate  al 
high  temperature  (80K  or  higher). 

The  counterdoping  technique  calls  for  doping  a  semiconductor  witf 
either  a  deep  donor  (Figure  1-1)  or  a  deep  acceptor  (Figure  1-2)  and  ther 
counterdoping  the  semiconductor  with  a  shallow  impurity  of  the  oppositr 
kind.  The  infrared  detector  thus  constructed  consists  of  an  infrarec 
active  center  (the  deep  impurity  [donor  or  acceptor])  which  has  los- 
charge  to  the  compensating  counterdopant.  The  deep  acceptor  is  nor 
negatively  charged  for  the  deep  acceptor  -  shallow  donor  case.  The  deei 
donor  is  positively  charged  if  we  have  the  deep  donor  -  shallow  accepto 
case . 

1  .2  BORON-DIVACANCY  COUNTERDOPED  SYSTEM 

We  have  chosen  the  shallow  acceptor  (0.045eV)  (Boron)  -  deep  dono 
(0.20eV,  0.73eV,  0.9eV)  (Divacancy)  system  to  construct  a  counterdope 
infrared  detector.  The  principal  infrared  active  dopant  is  the  dee 
donor.  The  concentration  of  boron  was  chosen  to  be  of  the  order  o 

3 

1  E  17  atoms/cm  ;  to  permit  the  deep  donor  to  be  the  active  element,  w 
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must  have  the  divacancy  content  be  of  the  same  order  of  concentration. 
It  is  desirable  to  have  the  number  of  divacancies  be  equal  to  or  slightly 
greater  than  the  number  of  boron  atoms. 

The  schematic  band  diagram  for  this  case  is  presented  in 
Figure  1-3.  We  assumed  for  this  case  that  +he  divacancy  is  either 
neutral  or  singly  positively  charged. 

"J  -  na  ■  na  for  nd  *  na  (1-1) 

Photoconductivity  occurs  when  we  have  photoabsorption  of  a  photon,  having 
energy  greater  then  Eg,  by  a  valence  electron.  The  donor  level  is  now 
neutralized  by  this  electron  and  a  hole  is  created  in  the  valence  band. 
The  holes  constitute  the  piiotocurrent ;  their  lifetimes  will  be  long 
because  the  recombination  capture  cross-section  of  holes  by  neutral 
donors  is  very  small.  Thus  the  operating  temperature  of  the  system  will 
be  higher  than  would  be  the  case  for  a  large  recombination  capture 
cross-section  detector. 

1 .3  PHASE  I  OBJECTIVES 

The  performance  of  our  Phase  I  activities  will  be  evaluated  against 
the  Phase  I  technical  objectives  provided  in  our  Phase  I  proposal.  We 
insert  Section  4.0  from  the  proposal  at  this  position.  When  we  discuss 
the  experimental  and  theoretical  results,  reference  will  be  made  back  to 
this  section. 


PHASE  I  TECHNICAL  OBJECTIVES 


An  experimental  and  theoretical  analysis  of  the  Bornn-Di vacancy 
counterdoped  infrared  detector  complex  will  be  undertaken  in  this 
Phase  I  research  effort. 

1.  Produce  an  identification  of  the  variables  which  influence 
the  BLIP  operating  temperature  for  the  Boron-  D-5 vacancy  counter 
doped  detector. 
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2.  Define  the  proper  annealing  procedures  and  temperature  to 
maintain  the  high  concentration  of  divacancies  produced  by  election 
beam  irradiation. 

3.  Characterize  the  levels  introduced  by  electron  beam 

irradiation.  Identify  crystalline  damage  defects. 

4.  Define  conditions  for  highest  photoconducti ve  response. 

5.  Develop  theoretical  model  for  the  rate  equations  to  predict 
proper  doping  concentrations  to  achieve  highest  BLIP  operating 
temperature . 

1  .4  DIVACANCY  PRODUCTION 

The  2MeV  electrons  supply  sufficient  energy  to  create  lattice 
vacancies  in  silicon  on  a  uniform  basis  in  depth.  The  threshold  energy 
for  displacement  of  a  silicon  atom  from  the  lattice  site  is  of  the  order 

n-2i 

of  30.0eV.  First,  four  bonds  are  broken  [E  =  4  times  (Bond 

Energy)  =  23.6eV].  Then  the  atom  goes  into  an  interstitial  position 

requiring  6.5eV  to  push  the  neighboring  atoms  out  of  its  path;  it  is 

necessary  to  surmount  the  barrier  set  up  by  these  atoms.  Finally,  all 

the  atoms  relax  about  the  vacancy.  Irradiation  was  performed  at  room 

temperature  where  the  lattice  vacancy  is  mobile.  This  provides  the 

longer  term  creation  of  divacancies.  It  is  necessary  to  permit  a  stable 

state  to  set  in  for  the  creation  of  divacancies  because  the  direct 

creation  of  divacancies  from  2MeV  electron  irradiation  is  only  about  55! 

n  -31 

that  of  the  creation  of  vacancies. 

The  rule  of  thumb  to  use  in  direct  creation  of  divacancies  fron 

-2 

electron  irradiation  is  to  state  that  10  divacancies  are  created  per 
electron . 
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2.0  EXPERIMENTAL  PROCEDURES 

2.1  FLOAT  ZONE  OPERATION:  BORON  DOPING 

A  Siemens  VZA3  vacuum  float  zoner  was  employed  to  purify  th 

polycrystall ine  silicon.  Four  zone  purification  passes  were  made.  The 

boron  was  added  on  a  fifth  pass  to  back  dope  to  a  boron  concentration  o 
3 

1  E  17  atoms/cm  . 

Each  zone  pass  of  the  coil  over  the  polycrystalline  rod  remove 
impurities  based  on  a  factor  termed  the  distribution  coef f icient( k) . 

C  =  Cq  [1  -  (1-k)  exp  (-kx/1)]  (2-1 

where  C  =  Concentration  of  element  in  zoned  rod 

Cq  =  Initial  concentration  of  impurity  element  before  pass 
*  =  Distribution  coefficient 

8.  =  Width  of  molten  zone 

x  =  Distance  from  seed  end  of  crystal 

Care  was  taken  in  removing  impurities  from  the  silicon  rod  becaus 
we  did  not  want  to  have  contribution  to  our  experimental  data  fro 
sources  other  than  boron  and  the  divacancy.  The  level  of  some  majo 
impurities  is  presented  in  Table  1-1  for  our  samples. 

TABLE  1-1.  IMPURITY  CONTENT  IN  SILICON 


Element 

Polysilicon  Concentration 
(atoms/cm3) 

Final  Single  Crystal 
Concentration  (atoms/c 

Boron 

2.6  E  12 

1  E  17 

Other  Group  III,  V 

5  E  12 

6  E  11 

Carbon 

5  E  16 

5  E  15 

Metals 

5  E  16 

1  E  13 
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Deep-Level  Transient  Spectroscopy  (DLTS),  which  utilize 

voltage  modulation  of  the  depletion  width  of  a  semiconductor  junction  ■ 
detect  electronic  transitions  at  defects  states  located  in  the  gap,  hi 
become,  during  the  past  ten  years,  the  most  widely  used  technique  fi 
microscopic  identification  of  electronic  defects  in  semiconductors.  Dl' 
exploits  nonradiative  transitions  between  localized  electronic  stati 
(deep  levels)  and  band  states  in  semiconductors,  and  it  thereto 
measures  thermal  transition  energies  which  may  determine  operation, 
behavior  of  semi  conductor  devices.  The  DLTS  is  the  most  sensiti' 
technique  to  detect  electrically  active  defects  and  provides  uniqi 
"signatures"  of  deep  levels.  DLTS  and  related  junction  spectroscop 
techniques  are  becoming  standard  methods  for  semiconductor  materia 
characterization,  processing  steps  monitoring  and  device  performam 
evaluation . 


The  essential  physics  of  DLTS  is  based  on  the  effects  of  trappii 

and  emission  of  carriers  at  deep  energy  levels  which  are  located  in  tl 

space-charge  layer  of  a  pn  junction  or  Schottky  barrier.  First,  \ 

consider  the  thermal-emission  transient  at  some  fixed  temperatun 
[2-2] 

Figure  2-1  shows  a  capacitance  transient  due  to  majority-carrii 

emission  with  insets  describing  the  condition  of  the  space-charge  1  ay< 
during  the  various  phases  of  the  transient.  A  capacitance  change 
caused  by  a  short  bias  pulse  which  introduces  carriers  and  thus  chang 
the  electron  occupation  of  a  deep  level  from  its  steady-state  value, 
this  population  returns  to  equilibrium,  the  capacitance  returns  to  i 
equilibrium  value.  If  we  now  consider  a  train  of  periodical  filli 
pulses  applied  on  the  reverse  biased  junction,  we  then  have  a  sign 
which  consists  of  a  series  of  transients  with  a  constant  repetiti 
rate.  As  the  temperature  is  varied,  the  time  constant  of  the  transien 
varies  exponentially  with  1/T.  This  is  shown  in  Figure  2-2  whi 
illustrates  the  basic  idea  of  the  DLTS,  the  so-called  rate-wind 
concept.  Thus,  the  DLTS  analysis  is  a  filtering  operation  and  there  a 
a  number  of  ways  in  which  the  DLTS  rate-window  concept  may  be  implement 
in  practice  and  a  DLTS  spectrum  produced.  The  first  method  proposed 
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[2-1  ] 

the  original  DLTS  experiments  by  D.V.  Lang  was  based  on  the  use 

a  double-boxcar  integrator.  The  DLTS  signal  is  the  difference  betwe 

the  transient  amplitude  at  two  times  t^  and  t^ .  after  the  pulse, 

shown  in  Figure  2-3.  There  is  zero  difference  at  these  two  gates  a 

consequently  no  DLTS  signal  for  either  low  or  high  temperatur 

correspondi ng  to  very  slow  or  very  fast  transients,  respectively.  T 

DLTS  peak  is  generated  when  the  transient-time  constant  t  is  on  t 

order  of  the  gate  separation.  The  rate  window  can  be  expressed  in  ten 

of  the  transient-time  constant  t  giving  rise  to  the  maximum  boxc 

"'H-ll 

output  as  a  function  of  temperature1  J,  namely. 


Tmax~ 


tl~t2 
1 n ( ti /t2 ) 


(2-2 


The  DLTS  scan  is  reversible  and  does  not  depend  on  the  magnitude  of  t 

heating  or  cooling  rate.  In  general,  the  DLTS  method  can  be  applied 

any  thermally  stimulated  process  which  allows  for  its  generation  by 

T  2-31 

pulsed  excitation.  J  The  DLTS  method  is  spectroscopic  in  the  sen 
that  the  various  DLTS  peak  heights  are  directly  proportional  to  t 
respective  concentration  of  individual  defects.  The  spectroscop 
character  of  DLTS  allows  for  quantitative  analysis  of  electrically  acti 

_5 

defects  with  a  sensitivity  of  the  order  10  of  the  shallow-lev 
doping  [NQ  -  Nft],  e.g.,  lO^cm  3  defects  if  |NQ  -  Nft|  =  1014cm  3.  The  DL 
technique,  when  applied  to  lightly  doped  material,  is  thus  the  mo 
sensitive  technique  to  detect  electrically  active  defects. 


The  activation  energy  of  a  deep  level  can  be  determined  in  t 
standard  manner  from  the  slope  of  an  Arrhenius  plot.  An  Arrhenius  pi 
is  constructed  from  data  obtained  by  recording  several  DLTS  spectra  usi 
different  rate  windows  and  by  plotting  the  log  of  the  rate  window  fr 
Reference  2-1  versus  the  inverse  temperature  of  the  DLTS  peak  T  \  1 

relationship  of  the  measured  activation  energy  to  the  true  depth  of  ' 
deep  level  is  complex.  This  is  due  to  the  temperature  dependence 
capture  cross-sections  and  due  to  the  temperature  dependence  of  • 
deep-level  energy  itself.^  3^ 
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Carrier  capture  measurements  can  be  done  by  recording  the  variation 

of  peak  height  as  a  function  of  pulse  width.  The  occurrence  of 

temperature-dependent  capture  cross-sections  must  be  taken  into  account 

[2-3] 

in  order  to  obtain  the  true  depth  of  the  level. 

2.3  THE  HALL  EFFECT 


The  concentration  and  Identity  of  electrically  active  Impurities  in 
silicon  may  be  measured  using  the  Hall  effect.  The  Hall  effect 

measurement  is  performed  by  passing  a  current  I  through  a  sample  which 

has  a  magnetic  field  H  applied  perpendicular  to  I.  The  Hall  voltage  V 

H 

is  measured  perpendicular  to  both  I  and  H. 

Initially,  the  magnetic  field  Is  turned  off;  for  H  =  0  we  find  that 

V  =  0.  When  the  magnetic  field  is  turned  on,  the  moving  charges, 
H 

represented  by  the  applied  current,  are  deflected  by  the  Lorentz  force  in 
a  direction  perpendicular  to  both  I  and  H.  At  steady  state,  the 

electrostatic  field  which  is  set  up  is  termed  the  Hall  field. 

The  Hall  coefficient,  R  ,  is  a  function  of  the  Hall  voltage,  the 

H 

current,  the  magnetic  field,  and  the  thickness  t  of  the  sample. 

RH  =  (VH/I)  (T/H)  (2-3) 

The  resistivity  of  the  sample  is  calculated  from  knowledge  of  the 
thickness  t,  the  width  w,  and  the  separation  of  the  voltage  contacts 
8. . 


P  =  (V/I)  (wt/2.) 


(2-4) 
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3 . 0  EXPERIMENTAL  DATA 


The  DLTS  spectra  exhibit  the  various  defects  which  are  created  by 

the  electron  irradiation.  We  see  in  Table  1-1  the  elements  of  interest 

which  may  be  involved  in  creating  defect  center  complexes.  The  boron 

3  . 

concentration  of  1  £  17/cm  is  so  high  that  it  is  likely  to  see  boron 
as  a  component  of  a  number  of  complexes.  A  concentration  of 

3 

5  E  15/cm  for  carbon  in  silicon  is  observed  as  a  best  value  for  any 
oil  pumped  float  zoned  silicon.  This  level  of  impurity  content  is 

sufficiently  high  to  produce  an  observable  center.  The  normally 
expected  vacancy,  interstitials  and  divacancy  defects  will  appear  in 
addition  to  (Boron-Vacancy) ,  boron  interstitial  and  carbon  interstitial. 

We  performed  the  electron  irradiation  at  300K;  the  vacancy  defect 
is  mobile  at  this  temperature  and  will  form  either  divacancies,  B  -V, 
or  other  complexes.  We  do  not  expect  the  DLTS  to  obcerve  the  vacancy 
defect  by  itself. 

3-1  DLTS  UNIT  NO.  1 

3.1.1  C-V  and  DLTS  Analysis 

Capacitance-Voltage  (C-V)  measurements  were  used  to  determine  the 

background  doping  and  Schottky  barrier  height  of  the  samples.  The 

barrier  height  measurement  gives  an  indication  of  the  quality  of  the 

diode.  On  p-type  material,  an  A1  Schottky  diode  would  usually  give  a 

barrier  height  of  about  1/3  Vgap  or  for  silicon,  0.34eV.  The  barrier 

2 

height  is  calculated  from  the  slope  of  the  1/C  vs.  V  curve  for  the 
diodes  fabricated  for  this  experiment  <$>  was  calculated  to  be 
0.16eV.  This  low  value  probably  indicated  a  reaction  of  the  A1  with  the 
surface  or  the  presence  of  surface  states. 

The  net  doping  concentration  calculated  from  the  C-V  data  was 
1 6  — 3 

3.8  x  10  cm  which  corresponds  well  with  the  value  of 

1 6  -3 

2.4  x  10  cm  calculated  from  Hall  effect  analysis.  C-V  data  are 
usually  higher  due  to  r-factor  corrections  to  the  Hall  data.  This  net 
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doping  concentration  is  used  to  calculate  the  trap  concentrations  in 
DLTS. 


The  C-V  data  were  taken  on  a  Hewlett-Packard  4280A,  1MHz,  C-V 

meter  under  control  of  an  LSI-11/23  computer  for  on-line  data 

2 

acquisition.  The  1/C  vs.  V  curve  is  plotted  in  Figure  3-1-1. 

3.1.2  DLTS  Equipment 

The  DLTS  system  is  a  fully  automated  system  consisting  of  an 
LSI-11/23  computer  with  a  30Mb  Winchester  disk  drive.  The  system  is 
shown  schematically  in  Figure  3-1-2.  The  AC  signal  is  derived  from  two, 
phase  locked,  frequency  synthesizers  which  allows  phase  shifting 
accuracy  of  0.1  degree.  Both  are  controlled  on  the  IEEE-488  bus.  The 
filling  pulse  is  provided  by  an  HP  1 1 6A  pulse  generator,  which  is  also 
IEEE-488  controllable.  The  "bridge"  used  to  balance  the  quiescent 
capacitance  was  fabricated  "in-house"  at  General  Electric  by  Dr.  Andrew 
Evwaraye. 

The  signal  is  processed  through  an  EG&G/PARC  115  preamplifier  and 
combined  with  the  phase  shifted  signal  from  the  HP3336C  at  the  mixer. 
( HP1 4541  A)  and  the  transient  signal  is  then  passed  through  an  EG&G15 
preamplifier  to  a  Tracor-Northern  TN-1710  signal  averages.  The  TH-1710 
captures  the  transient  and  outputs  rate-window  data  to  the  LSI-11/23. 
Up  to  twenty-five  rate  windows  may  be  programmed  simultaneously  on  the 
TN-1710.  The  data  is  analyzed  in  a  real-time,  interactive  program  on  a 
PRIME  850  computer. 

Data  may  be  acquired  at  any  frequency  between  lOKHz  and  20MHz. 
The  data  discussed  in  these  results  was  taken  at  20MHz. 

3.1.3  Experimental  Results 

Three  peaks  were  detected  in  the  sample  examined.  These  peak 

positions  are  associated  with  a  t  =1.9  E-3.  The  activation 

max 


14 


IJES00I  .  CVDRTH 


15 


•B.sa  H.  19  I.M  l.Sfi  <? .  2S  ?.9»  J.C3  1 .  J 1  b.W 

VDL  ffi'A  (V) 


OUTPUT 


Ficure  3-1-2  Eauipment  Block  Diaqran 


15 


energies  calculated  were  0.1  geV,  0.35  eV,  and  0.54  eV  respectively. 
The  acccepted  literature  value  for  the  divacancy  is  0.21  ±  .02eV. 

The  activation  energies  are  calculated  from  an  Arrhenius  plot 

2 

Fiqure  3-1-3  of  In  (T  t  )  vs  1/KT  by  a  linear  regression 

3  max  — 

fit.  The  spectrum  is  shown  in  Figure  3-1-4.  The  shoulder  on  the  left 
of  the  V-O-C  is  probably  due  to  the  b-C  complex,  but  could  not  be 
resolved  well  enough  to  calculate  an  energy.  This  peak  is  diminished 
due  to  the  annealing  of  the  sample  after  irradiation.  The  Arrhenius 
plot  is  shown  in  Figure  3-1-4. 

3.2  DLTS  UNIT  NO.  2 

Defects  in  2MeV  electron-irradiated  silicon  samples  were  analyzed 
using  a  Polaron  DL4600  DLTS  system  over  a  temperature  range  of  90-300K. 
All  DLTS  scans  were  performed  in  darkness  with  decreasing  trmperature. 
The  samples  were  analyzed  prior  to  electron  bombardment  to  determine  the 


background 

residual  impurities. 

No 

defect 

levels  were 

found 

at 

12 

-3 

then 

concentrati ons 

exceeding  10 

cm  . 

The 

samples 

were 

i rradiated 

with 

2MeV  electrons 

near 

room 

temperature 

( 300K) 

to 

2xl017/cm2. 

The  thermal  stability  of  the  defects  was  determined  by  isochronal 
annealing,  with  10-minute  intervals,  up  to  450K  in  situ  under  a  vacuum 
of  10  2  torr.  Dry  nitrogen  atmosphere  was  used  for  higher 
temperatures ,  and  the  DLTS  spectrum  was  measured  after  each  stage  of 
anneal ing. 

We  observed  three  distinct  defects  of  comparable  concentration 
that  are  produced  by  room-temperature  electron  irradiation  in 
boron-doped  FZ  silicon.  We  tentatively  correlated  the  high 
temperature  trap  at  Ey  +  0.44eV  with  the  carbon-interstitial  + 

v  r  3-1 1 

carbon-substitutional  (C.  +  C^)  pair  identified  by  Brower, 
using  EPR  techniques,  in  e  1  ectron-'- ^radiated  boron-doped  silicon  which 
contains  carbon.  Brower  suggests  that  the  pair  is  formed  when  a  mobile 
carbon  interstitial  is  trapped  by  a  substitutional  carbon  atom.  It  can 
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be  produced  by  either  electron  or  neutron  irradiation  at  room 

[3-21 

temperature.  The  pair  was  observed  by  Jaworowski  et  al.  at 

Ev  +  0.39eV  as  the  major  DLTS  peak  after  electron-i rradiation  of 
silicon  with  high  concentration  of  substitutional  carbon.  Our 
correlation  is  based  on  the  similar  behavior  under  isochronal  annealing 
showing  the  center  disappearing  at  approximately  575K. 

Figure  3-2-1  shows  typical  DLTS  spectra  for  the  two  low 

temperature  defects  H  and  measured  with  six  different  rate 

windows.  Their  thermal  activation  energies  determined  in  the 

+ 

conventional  manner  (Figure  3-2-2)  were  Ev  +  0.15  0.02  eV  and  Ey  h 

0.28  ±  0.02  eV  and  they  were  detected  in  concentrations  higher  thar 
15-3 

2  x  10  cm  .  Capture  cross-sections  were  not  determined,  anc 
additional  possible  corrections  for  their  temperature  dependence  are  nol 
included.  Figure  3-2-3  shows  results  of  isochronal  annealing  of  anc 
traps.  Both  traps  anneal  out  rapidly  between  500  and  550K. 

The  tentative  identifications  given  here  are  based  on  comparisoi 
with  energy  levels  and  the  thermal  stability  of  defects  previously 
observed  in  boron-doped  silicon  which  had  been  electron-i rradiated  a> 
near  room  temperature.  The  H  deep  level  corresponds  to  a  positive 

+  .  r  3-3 1 

state  of  the  divacancy  (V^).  KimerlmgL  J  first  correlated  the 
H  peak  with  a  positive  charge  of  the  divacancy  which  was  firmly 

r  3—4 1 

identified  by  Watkins  and  Corbett1-  J  in  their  electron  paramagnetii 

resonance  experiment.  The  defect  is  most  likely  a  boron  relatei 
trap.  If  the  impurity  concentration  is  high,  the  mutual  anihilation  o 
the  primary  radiation  defects,  i.e.,  vacancies  and  interstitials 
decreases  and  the  formation  of  complexes  with  the  basic  dopan 

increases.  The  activation  energy  of  the  H  is  equal  to  th 

r  3—51  2 

vacancy-oxygen-boron  complex  ( V— 0-B ) , L  J  which  can  be  observed  in  F 
silicon  after  illumination  at  room  temperature,  but  the  (V-O-B)  comple 
anneals  out  at  675K  which  is  approximately  125K  higher  temperature  tha 
the  annealing  temperature  of  the  defect.  Chantre1  J  recent! 

suggested  the  possibility  of  a  room  temperature  stable  boron-vacanc 
(B-V)  pair.  This  is  a  new  configurationally  bistable  center  in  quenche 
boror,  doped  silicon.  The  observed  thermal  stability  of  the  center  (day 
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at  room  temperature)  seems  to  be  in  conflict  with  the  observed  annealing 

temperature  for  the  (anneals  out  above  500K).  The  B-V  pair  may  be 

unstable  at  room  temperature  in  the  presence  of  the  high  density  of 

electrons  and  holes  generated  by  2MeV  electron-irradiation.  This  would 

explain  the  nondetection  of  the  bistable  center1  J  in  room 

temperature  electron-i rradiated  boron-doped  silicon.  A  prominent 

f  3-71 

majority  earner  tn.p  was  observed  by  Drcvinsky  and  CcAnyelii 

after  electron-electron  irradiation  in  FZ  silicon  heavily  doped  with 

1 7  -3 

boron  (10  cm  )  at  Ey  +  0.30eV.  They  suggest  a  boron 
interstitial-boron  substitutional  pair  (B  -B  )  as  a  candidate 

f3_31  s 

structure.  J  The  annealing  behavior  of  their  center,  however,  is 

r  3-5] 

similar  to  the  (V-O-B)  center;  t . . e i r  center  does  not  anneal  out 

until  675K . 

1 7  -3 

The  doping  impurity  level  (B:  2  x  10  cm  )  in  our  FZ  silicon 

15  -3 

is  much  higher  than  the  oxygen  concentration  (0:  2  x  10  cm  ). 


This  supports 


possible  assignment  of 


to  the  boron 


interstitial-boron  substitutional  (B.-B  ).  A  new  deep  level  has 

is 

appeared  after  the  5 00K  annealing  step  and  annealed  out  at  550K 

(Figure  3-2-3).  Similar  deep  levels  were  observed  in  boron-doped  silicon 

r  3—5 1 

after  the  annealing  of  the  divacancy  at  about  550K. 

r 3 _ 8  3-91 

Recently,  Londos  and  Banbury1  ’  J  reported  on  DLTS  and  MCTS 
(minority-carrier-trap  spectroscopy)  studies  for  electron-irradiated 
boron-doped  silicon.  Song  et  al.^3  observed  using  DLTS  a 

configurational  metastability  of  vacancy-related  defects  in 
electron-irradiated  n-type  silicon.  These  recent  studies  have  raised  a 
number  of  important  new  questions  concerning  the  properties  of  primary 
radiation  defects  and  their  interactions  in  silicon. 


3.3  HALL  ANALYSIS 


A  1cm  by  1cm  sample  was  taken  from  the  electron  irradiated  silicon 

3 

sample.  The  boron  concentration  was  1  E  17  atoms/cm  .  Divacancies 
were  introduced  by  electron  bombardment  (2MeV).  We  assumed  a  divacancy 


24 


creation  factor  of  0.01  divacancies  for  every  electron.  The  intent  was 
to  produce  as  many  di vacancies  as  we  had  boron  atoms. 

Wires  were  indium  soldered  to  the  four  corners  of  the  1cm  by  1cm 
sample.  This  sample  was  intended  to  be  used  later  for  low  temperature 
infrared  absorption  measurement;  thus,  no  laser  annealing  was  performed 
on  the  four  contacts. 

Hall  data  were  taken  from  23K  to  380K.  The  low  temperature 
results  demonstrated  some  scatter  in  the  data  and  poor  r^/r^  values 
because  we  did  not  laser  anneal  the  contacts.  Because  of  this  problem, 
we  do  not  include  data  between  23K  and  26K.  If  possible,  at  a  later 
date,  we  shall  redo  this  temperature  range  with  laser  annealed  contacts. 

The  sample  was  p-type  over  the  entire  temperature  range; 
indicating  that  the  number  of  divacancies  created  did  not  equal  the 
number  of  boron  atoms  present.  The  hole  concentration  as  a  function  of 
1 000/T  is  presented  in  Figure  3-3-1. 

The  data  were  analyzed  by  fitting  the  charge  balance  equation. 
The  high  temperature  data  above  200K,  which  indicate  the  presence  of  a 
deeper  level,  were  deleted  for  the  initial  fits.  From  these  fits,  we 
conclude  the  following  concentrations  for  the  elements  of  interest. 

3 

Donor  Concentration  =  7.36  E  15/cm 

3 

Boron  Concentration  =  2.40  E  16/cm 

3 

2nd  Acceptor  Level  Concentration  =  9.61  E  15/cm 
The  activation  energy  of  this  second  acceptor  level  is  0.034eV. 

The  calculated  fit  using  these  values  is  presented  in  Figure  3-3-2. 

The  next  set  of  calculations  used  all  the  data  to  provide  a  fit. 
The  parameters  for  donors  and  the  first  two  acceptor  levels  did  not 
change  significantly.  An  activation  energy  of  0.182eV  was  indicated  for 
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Figure  3-3-2  Calculated  Fit 


the  deep  level.  This  is  within  experimental  error  of  the  energy  of 

0.21eV  generally  associated  with  the  divacancy. 

We  performed  another  set  of  calculations  by  fixing  the  energy  of 

2 

the  divacancy  at  0.21eV.  In  this  case,  the  fit  produced  a  x  ,  which 

is  a  measure  of  how  well  the  data  have  been  fit,  of  0.84.  This  compares 
2 

with  a  x  of  0.75  when  all  the  parameters  were  allowed  to  vary.  Both 
these  fits  may  be  considered  to  be  equivalent. 

The  last  fit  which  set  the  divacancy  level  at  0.21eV  produced  the 
trace  shown  in  Figure  3-3-3.  The  parameters  associated  with  this  trace 
are  listed  below. 


N 

N 

E 

N 

E 

N 

E 


D 

B 

B 

A2 

A2 

A3 

A3 


=  6.96  E  15/cm3 
=  2.47  E  16/cm3 
=  0.0469eV 
=  8.98  E  15/cm3 
=  0.0339eV 
=  3.14  E  16/cm3 
=  0.21eV  (Fixed) 


3.4  INFRARED  ANALYSIS 


The  infrared  spectroscopic  technique  permits  examination  of  the 
defects  produced  in  silicon  by  2MeV  electron  irradiation.  Interstitial 
and  vacancy  defects  extend  the  bandedge  into  the  infrared  when  the  sample 
is  observed  at  300K.  These  defects  cause  the  normally  sharp  1.06  micron 
bandedge  to  tail  into  the  infrared  to  1.4  microns.  This  tail  can  be 
annealed  out  only  at  1200K  heat  treatment. 

Observation  at  300K  also  produces  a  band  at  1.8  microns  due  to  the 
neutral  divacancy  and/or  the  negative  divacancy.  When  the  Fermi 
level  is  below  Ec  -  0.22eV,  the  negative  divacancy  contributes 
to  the  1.8  micron  band. 
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~e  3-3-3  Calculated  Fit  =  0.21eV 


When  the  Fermi  level  is  above  Ec  -  0.22eV,  the  divacancy  state 

shifts  over  to  the  double  negative  charged  state  and  may  be 

observed  at  3.45  and  3.62  microns.  Temperatures  below  60K  must  be  used 

to  observe  this  state;  the  bands  become  sharper  and  more  intense  as  the 

temperature  is  decreased  below  60K.  The  position  of  the  Fermi  level  will 

determine  the  relative  strength  of  the  1.8  micron  and  the  3.45,  3.62 

micron  bands.  As  stated,  the  Fermi  level  must  be  above  E  -  0.22eV  to 

c 

populate  the  double  negative  charged  state.  Thus,  to  observe  the 
strongest  photoconductivity  signal  at  the  4  micron  band  position,  this 
condition  must  be  met. 

The  infrared  analysis  presented  for  our  electron  irradiated  samples 
can  not  provide  an  exact  value  for  the  boron  concentration  because  the 

3 

concentration  is  too  high  (1  E  17  atoms/cm  ).  The  far  infrared  curve 
of  Figure  3-4-1  shows  the  reversal  of  the  spectral  lines  indicating  the 
high  boron  content.  Phosphorus  content  can  normally  be  observed  by 

illuminating  the  sample  in  the  course  the  scan.  No  phosphorus  was 

observed;  this  indicates  that  the  multipass  zoning  did  successfully 

remove  phosphorus.  There  can  be  no  lower  level  for  the  phosphorus 

content  because  the  boron  content  was  so  high.  Normally  we  desire  low 
levels  for  boron  to  obtain  a  good  limit  of  detection  statement  of 
phosphorus . 

The  main  curve  of  interest  is  produced  in  Figure  3-4-2.  This 

presents  evidence  of  the  V*  center  for  the  divacancy  at  2500 

_1  £- 

cm  .  This  spectra  also  was  produced  at  a  temperature  of  5K.  The 

concentration  of  divacancies  was  not  greater  than  the  boron 

concentration;  however  we  were  successful  in  observing  the  defect 

structure  desired.  The  divacancies  which  were  produced  were  compensated 

by  the  boron  in  the  crystal  and  thus  exposed  the  \L  state.  The 

-3  ^ 

boron  concentration  was  1  E  17  cm  which  was  about  25  times  the 

divacancy  concentration.  The  DLTS  analysis  gave  a  boron  value  of 

-3 

3.8  E  16/cm  for  boron  and  the  Hall  data  gave  a  value  of 

-3  -3 

2.4  E  16/cm  .  This  compares  with  an  infrared  value  of  1  E  15/cm 

for  the  state. 
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The  vj  state  observed  is  the  desired  level;  it  results  from 
the  compensation  of  the  boron  atoms.  We  do  not  observe  neutral  or 
negatively  charged  divacancies  because  boron  has  performed  its  desired 
role  of  compensation.  The  intensity  of  the  positive  divacancy  would  be 
greater  if  we  had  a  great  dose  of  electrons.  The  desired  level  would  be 
an  equal  concentration  of  divacancies  and  boron  atoms.  The  interesting 
fine  tuning  experiments  in  which  the  concentrations  are  varies  about  this 
exact  compensation  value  are  beyond  the  scope  of  this  Phase  I  effort. 

An  a+tempt  was  made  to  illuminate  the  sample  at  the  same  time  the 
spectra  was  taken.  There  was  no  change  in  the  intensity  of  the  positive 
divacancy  and  no  neutral  divacancies  or  negative  divacancies  were 
created.  This  provides  important  evidence  that  Lhe  counterdoping  theory 
of  Elliott^  ^  is  correct;  infrared  detectors  may  be  created  by 
counterdoping  and  are  stable  under  illumination. 


4.0  APPLICATION  OF  MULTIVALENT  STATISTICS  TO  COUNTERDOPED  DETECTORS 


4.1  INTRODUCTION 

Host  impurities  that  are  contemplated  for  their  use  as  either  deep 
donor  or  acceptor  in  the  counterdoped  infrared  detector  supply  several 
other  energy  levels  in  the  gap.  The  multivalent  nature  of  these  centers 
has  never  been  incorporated  in  a  study  of  the  expected  detector 
response.  In  particular,  these  levels  may  be  ordered  in  the  usual 
donor-acceptor  ascending  order  of  energy  levels;  or  they  may  be  reverse 
ordered  as  is  character! stic  of  the  negative  -U  center.  The  consequences 
of  the  level  ordering  and  the  multivalent  nature  of  these  impurities  is 
the  subject  of  the  following  theoretical  analysis. 

4.2  THE  DIVACANCY 

The  divacancy  in  silicon  introduces  several  levels  in  the  energy 
gap  of  the  host  lattice.  In  p-type  silicon  the  first  level  nearest  the 
valence  band  is  the  donor  level  associated  with  the  charge  transition. 

+  0 

V2  -  V2  ,  (4-1) 

which  may  occur  with  the  emission  of  a  hole  to  the  valence  band,  or  by 
capture  of  an  electron  from  the  conduction  band.  The  second  level  of 
interest  is  the  level  corresponding  to  the  transition. 

0 

V2  -  V2  .  (4-2) 

The  divacancy  can  also  carry  one  more  negative  charge  in  the  transition. 

V2  -  V2  .  (4-3) 

All  these  charge  states  are  realizable  in  equilibrium  as  the  Fermi  level 
is  raised,  or  in  a  non-equilibrium  experiment  via  carrier  injection. 

Figure  4-1  illustrates  the  present  state  of  knowledge  about  the  divacancy 

,  ,  .  .,.  [4-1] 

levels  in  silicon . 
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The  transition  in  Equation  4-1  is  important  for  the  operation  of  the 

counterdoped  detector.  It  is  of  interest  to  study  the  effect  of  the  next 

charge  state,  Vp  in  Equation  4-2,  on  the  detector  performance.  The 

model  used  so  far  has  included  only  the  lowest  donor  level  of  the 
f  4-21 

divacancy.  Here,  the  next  divacancy  level  will  be  considered  and 

incorporated  into  the  model.  For  the  sake  of  generality,  the  initial 
developments  will  apply  to  any  counterdoped  system  with  two  charge  states 
in  the  gap.  Reference  4-2  is  included  as  Appendix  1  for  your 
information.  The  theory  developed  in  this  reference  has  a  bearing  on  the 
total  picture  of  counterdoped  detectors.  Extended  experimental  research 
is  necessary  to  test  the  predictions  presented  in  this  theory. 

4.3  MULTIVALENT  STATISTICS 


We  must  start  first  with  the  statistics  of  multivalent  levels. 

F4-31 

This  subject  has  been  covered  extensively  by  Shockley  and  Last, 

,  [4-4]  [4-5] 

Blakemore,  and  Look.  Let  the  state  of  the  impurity  with  no 

ionizable  electrons  be  the  zeroth  state  of  the  impurity.  This  state  will 

have  the  total  energy  which  will  be  set  to  zero,  cQ  =  0,  and  the 

degeneracy  yo.  The  state  with  k  ionizable  electrons  will  be  denoted 

with  the  subscript  k;  and  y^,  for  the  total  energies  and 

degeneracies,  respectively.  In  the  case  of  the  divacancy  in  silicon,  the 

k  =  0  state  is  the  positively  charged  divacancy  V*  the  k  =  1  state 

0  d 

is  the  neutral  divacancy  V^,  etc.  The  energy  needed  to  change  the 

charge  state  +/0  is  -  c^.  This  energy  will  promote  an 

electron  from  the  valence  band  to  the  first  level  in  the  gap  (the  0.204eV 
level  in  Figure  4-1 ) . 

It  is  understood  that  the  system  cannot  support  any  other  charge 

states  in  the  gap.  For  example,  the  triply  negative  charge  state 

V~  would  be  a  resonance  in  the  conduction  band  and  as  such  an 

autoionizing  system.  The  extra  electron  would  be  at  the  bottom  of  the 

conduction  band  rather  than  on  the  divacancy.  Similarly,  the  double 

positively  charged  state  V^+  is  a  resonance  in  the  valence  band,  or 
below,  so  that  it  takes  more  than  the  bandgap  energy  to  remove  an 

electron  from  the  V*  charge  state.  That  is  why  V*  has  no 
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ionizable  electrons  and  why  it  takes  no  energy  to  transfer  one  electron 

from  the  valence  band  to  form  v*.  By  the  same  token,  simple  donors 

in  silicon  have  either  four  or  five  outer  shell  electrons  in  the  D  and 

charge  states  respectively.  Only  one  electron  is  ionizable  though 
0 

in  the  D  configuration. 

The  occupancy  of  N  impurities  with  k  ionizable  electrons  can  be 
derived  most  simply  by  the  maximum  probability  method.  The  number  of 

impurities  with  k  ionizable  electrons  is  given  by  Equation  B26  of 

,  [4-5] 

reference. 


N 

nK  =  -  (4-4) 

1+s  (yl/yk)  exp  [ck-cj,-(k~i)v]/kT 


where  =  chemical  potential 

k  =  Boltzman  constant 

T  =  absolute  temperature 


As  emphasized 

by  Look^4  5^,  the 

relevant  parameters 

are 

the 

relative  degeneracy 

factors,  y  /y.  , 

and  the  differences 

in 

total  energies, 

-  e  .  The  relevant  parameters  are 

not 

the 

single  electron  energy 

level s . 

Equivalently,  the  ratio  of  concentrations 


nk  T*  ( v-ck ) 

-  =  -  exp  - 

nk— 1  Tk-1  kT 


( 4— 4a ) 


where  c  is  the  energy  required  to  take  the  *cth  electron  from  the 

1C 

valence  band  and  place  it  onto  the  defect  with  k-1  electrons  already 
present,  i.e., 


Ek  =  ck  -  ck-l  =  E ( *c/k-1  ) 


( 4  —4  b ) 


The  new  symbol  E(«/tc-l)  gives  the  "before"  and  "after"  picture  of  the 
number  of  ionizable  electrons  on  the  defect.  We  shall  also  use  the  symbol 
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9k  =  Yk/Yk-1 


(4-4c ) 


for  the  "relative"  degeneracy  of  a  level. 

Exploration  of  the  detector  response  of  a  multivalent  impurity  system 
offers  an  intriguing  possibility  of  contemplating  the  behavior  of  a 
negative  -U  system.''4  &  to  4  ^ese  systems,  the  total  energy  of 

the  k+1  electron  configuration  is  lower  than  the  configuration  with  k 
electrons . 

ck+l  <  ck  •  (4~5) 

It  is  as  if  the  extra  electron  was  attracted  rather  than  repelled  from 
the  k  electrons  already  on  the  impurity.  In  fact,  it  is  the  Jahn-Teller 
relaxation  of  the  defect  in  the  (k+1)  charge  state  which  lowers  the 
electronic  energy  while  raising  the  elastic  energy.  The  overall  gain  in 
energy  is  negative,  hence  the  negative  -U  property  of  the  center. 

4.4  CHARGE  BALANCE  FOR  THE  DIVACANCV 


In  the  case  of  the  divacancy  system,  the  charge  balance  equation 
will  be 


where 


n  +  N'  +  2N=  +  =  p  +  Nfr  +  N+ 


(4-6) 


N 

N' 

N= 

P 

N. 


densi 

n  = 
o 

nl  = 

n2  = 

n3  = 
densi 

densi 

densi 


ty  of 
densi 
densi 
densi 
densi 
ty  of 
ty  of 
ty  of 


electron  in  the  conduction  band 
ty  of  positively  charged  divacanci 
ty  of  neutral  divacancies 
ty  of  negatively  charged  divacanci 
ty  of  doubly  negative  charged  diva 
holes  in  the  valence  band 
the  counterdopant,  here  boron 
the  shallow  donors 


es 

es 

cancies 
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The  conservation  of  divacancies  accounts  for  the  total  density  of 
di vacancies 


N  =  N+  +  N°  +  N"  +  N" 


(4-7) 


4.5  DEGENERACY 


Next,  we  must  determine  the  degeneracy  of  the  divacancy  charge 
states  and  their  total  energy  for  input  into  Equations  4-4,  4-6,  and 
4-7.  The  easiest  task  is  to  determine  the  relevant  total  energies  since 
these  are  experimentally  measured.  In  the  notation  of  Equation  4-4  and 
Figure  4-1 , 


g  =  0  ,  (4-8) 

o 

^  =  E1  =  0.204eV  =  Ep  , 

*2  -  E,*2  -  0.93ev  -  Ed+Eda  . 

c3  =  ET+E2+E3  =  1  R°V  =  Ed+Eda+Edaa  , 
where  E.  is  the  energy  necessary  to  add  the  *th  electron  from  the 

IN 

valence  band  to  the  (k-1)  charge  state  of  the  divacancy. 

Explicity  we  have  using  Equation  4-4, 

(4-9a) 

(4-9b) 

( 4—9  c ) 


N 


no  = 


Y1  (u  -  El)  Y2  r(2v  -  El  -  E2) 

1+ -  exp  [ - ]+ - exp[ - ] 

YO  YO  kT 


=  N+ 


N 


ni  = 


.  YO  r(El-p),  r  (  2 )  , 

1  + exp[ — ]+ exp[ — — — ] 

Y1  kT  Y1 


=  nO  , 


N 


no  = 


YO  r  (  E t  -hE 2_ 2 vi )  yi  (E2-p) 

1+ _ exp[ - — - ]+ - exp[ - —] 

Y2  kT  Y2  kT 


=  N"  , 
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(--/-)  O.c 

(-/O)  0.7; 

(0/+)  0.2 


Divacancy  in  Silicon 


where  the  N  =  n3  state  has  been  dropped  out  for  physical  reasons 

(low  occupancy  level  at  the  temperatures  of  interest). 

The  degeneracy  factors  in  Equations  4-4c  and  4-9,a,b,c  depend  on 
the  detailed  electronic  structure  of  the  divacancy.  Imagine  two  nearest 
neighbor  silicon  atoms  in  the  diamond  structure.  Eacn  of  the  atoms  have 
three  other  silicon  nearest  neighbor  atoms.  Bonds  connecting  every  pair 
of  atoms  have  two  electrons  since  each  atom  contributes  four  electrons 
for  bonding.  When  two  atoms  are  removed,  a  divacancy  is  created.  The 
bonds  with  six  remaining  atoms  are  broken,  and  there  are  six  dangling 

bonds  pointing  toward  the  position  of  the  divacancy.  In  the 
Watkins-Corbett  model, ^  ^  we  find  two  pairs  of  atoms  rebonding;  this 
ties  up  four  electrons  and  leaves  two  loosely  bound  electrons 
representing  the  neutral  state  of  the  divacancy,  V^.  The  bonding 
state  of  these  loosely  bound  electrons  accommodates  either  one  electron 

in  the  \L  charge  state,  or,  as  mentioned  above,  two  electrons  in  the 
0  ^ 

state.  Two  additional  electrons  can  be  placed  in  the  antibonding 
combination  of  the  two  dangling  orbitals  to  produce  in  succession  the 

and  V?  charge  states. 

The  symmetry  point  group  of  the  undistorted  divacancy  would  be 

0o,.  The  Jahn-Teller  distortion  lifts  the  electronic  degeneracy, 

3d 

resulting  in  a  rebonding  of  the  two  atom  pairs;  thus,  the  symmetry  is 
reduced  to  that  of  the  C  ^  point  group.  The  degeneracy  of  each  level 

is  two-fold  counting  the  spin.  The  energy  level  diagram  for  the  four 
charge  states  is  given  in  Figure  4-2. 

In  Figure  4-2,  the  number  of  ways  the  given  number  of  electrons  can 
be  placed  gives  the  degeneracy  factor.  In  order  to  avoid  confusion,  note 

4“ 

that  the  lone  electron  in  the  state  is  not  an  ionizable 

electron.  It  takes  more  energy  to  ionize  it  than  to  promote  a  valence 

electron  to  the  conduction  band. 

The  V*  and  states  have  unpaired  electrons  and  are 

Electron  Paramagnetic  Resonance  (EPR)  active.  They  were  observed  by 

Watkins  and  Corbett.^  There  are  newer  models  of  the 
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Figure  4-2  Electron  Occupancy  of  the  Dangling 
Bond  States  of  the  Divacancy  and 
the  Divacancy  and  the  Associated 
Degeneracy  Factors 
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divacancy;^4  ^  however,  much  work  remains  to  elucidate  all  of  its 
properties.  The  Watkins-Corbett  model  is  sufficient  for  the  accuracy  of 
the  present  calculations. 

We  shall  define  the  degeneracy  factor  g  as 


For  simple 
For  simple 

4 . 6  RATE  EQUATIONS 


9D  = 

T1 

=  1/2  , 

TO 

9DA  = 

T2 

=  1 

TO 

9DAA  = 

Y3 

=  1/2  , 

TO 

acceptors 

TO  = 

4,  Yl  =  1 

donors 

TO  = 

1.  Tl  =  2 

(4-10a) 

,  ( 4-1  Ob) 

(4-10c) 

thus  g  =  1/4 
thus  g  =  2 


It  will  be  assumed  that  the  first  three  charge  states  of  the 
impurity  will  be  sufficiently  occupied  at  the  temperature  of  interest 
(near  liquid  nitrogen  temperature).  This  also  implies  that  the  two 
electronic  levels  involved,  0/+  and  -/0,  can  communicate  only  with  the 
valence  band.  The  rate  at  which  the  density  of  the  charged  defects 
changes  is  given  by 


—  =  B(0/+)N0  p-A(0/+)N+  (4-11) 
dt 

The  first  term  is  due  to  hole  capture  by  and  the  second  term  is  the 
result  of  hole  emission  to  the  valence  band.  The  B  coefficient  is  the 
recombination  coefficient  for  the  indicated  charge  state  transformation. 
The  A  coefficient  is  the  generation  coefficient  for  neutral  impurities. 


=  A(°/+)N+-B(°/+)N°  p-A(-/0)N°+B(~/0)N-  p  .  (4-12) 

at 
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The  first  two  terms  are  as  defined  in  Equation  4-11.  The  third 
term  is  the  result  of  hole  emission,  while  the  fourth  term  is  due  to  hole 
capture  for  the  negatively  charged  species. 


3N~ 

at 


=  A(  /0)N^-B(-/0)N~  p+(neglected  terms)  . 


(4-13) 


The  neglected  terms  involve  (“/-)  conversion. 

These  equations  are  supplemented  by  the  charge  neutrality  condition, 

p  +  N+  =  N§  +  N“  ,  (4-14) 

where  N§  is  the  density  of  a  shallow  counter  dopant,  in  this  case 
boron.  The  "conservation"  of  divacancies  completes  the  picture. 

N=N++N°+N"  .  (4-15) 

The  time  derivatives  in  Equations  4-11  to  4-13  are  zero  for  the 
steady  state  condition;  thus,  we  may  use  Equations  4-11  to  4-13  to  obtain. 


A0  =  A(0/+)=B(°/+)  —  P  ,  (4-16) 

N+ 

ADA  E  A(-/0)=B(-/0)  —  P  .  (4-17) 

NO 

In  thermal  equilibrium,  we  have 


P  =  NV  exp  (-v/kT)  (4-18) 

where  Ny  is  the  effective  density  of  states.  Using  Equation  4-4a, 
we  see  that 


N°/N+  =  n-|/n0  = 

(T°/Y+)exp[  — 

‘  EDA)n 

J  * 
kT 

(4-19) 

N~/N^  =  n2/n^  = 

(yT/y°)  exp 

v-^DA 

(4-20) 

kT 
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Thus  the  relationship  between  the  thermal  generation 

and  the  recombination  coefficient  B  become. 

coefficients  A 

-En 

Aq  =  Ny  Bp  ggexp( - ) 

kT 

(4-21) 

Eqa 

aDA  =  NV  bDA  9DA  exP  ( - )  • 

kT 

(4-22) 

Using  Equations  4-21  and  4-22  we  may  now  rewrite  Equations 

4-11  to  4-13  and  add  the  optical  generation  terms.  Equation  4-11  picks  up 

-  <oq  4>>  Nft  , 

Equation  4-12  picks  up 

(4-23) 

<oQ  4>>  N$-<o0a  4>>  Nq  0 

Equation  4-13  picks  up 

(4-24) 

<oqa  4>>  Ng 

(4-25) 

The  o'  s  are  photon  absorption  cross-sections  and  <t> 

flux.  The  brackets  are  spectral  integrals. 

It  is  convenient  to  define 

is  the  photon 

aD=gD  Nvexp(-ED/kT)+<oD  4»/BD  , 

(4-26) 

QDA=9DA  Nvexp(-E0A/kT)+<°0A4)>/BDA 

(4-27) 

In  terms  of  quantities  just  defined  the  steady-state 

of  the  charge  states  of  the  divacancy  are 

concentrations 

N*  -  N  f2 

p2-Hap  p+aD  aDA 

(4-28) 

N°  -  a°  P  N 

p2+aD  P+aQ  aDA 

(4-29) 

aD  aDA  N 

N  =  - 

P^+aoP+a0  aDA 
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(4-30) 

In  fact,  these  equations  are  simply  related  to  Equation1.  ,-9a  to  4-9c. 
The  solution  of  Equations  4-28  to  4-30  depends  on  the  density  of  holes  p 
at  the  given  temperature.  That  comes  from  the  charge  balance 
Equation  4-14  supplemented  by  the  equation  for  the  density  of  ionized 
acceptors  (Fquatio;.  4-4a), 


Na/n8  =  (YA/T^exPL— - — ] 


(4-31) 


kT 


The  acceptor  level  is  two-folo  degenerate  with  three  electrons  in 

it,  when  it  is  neutral.  It  is  nondegenerate  when  it  has  four  electrons 

0 

and  is  negatively  charged  (ya  =  1 , y A  =  4).  Thus 


NA  = 


1+4  exp [  (  EA  -  y)  AT] 


A  rate  equation  for  the  counterdopant  gives  rise  to 

&A  nA 


Na  = 


P  + 


where 


(Nv) 


exp  -  (EAAT)  +  <  <?a4>>/Ba  . 


(4-32) 


(4-33) 


Using  the  charge  balance  equation  (Equation  4-17),  we  obtain  a 
quartic  equation  for  p 

P4  +  p3  (a a  +  aD  +  N)  +  P2  [°D  (“A  +  aDA)  +  aA  (N  "  na)1  +  (4-34) 

p  aD  (aDA  aA  -  N  aDA  “  NA  aA)  "  aDA  aD  aA  (NA  +  N)  =  0 

If  we  set  a..  =  0,  we  obtain  the  equation  for  the  case  of  a  singly 

UA  14-21 

ionizable  deep  donor.  J 


4.7  DETECTOR  RESPONSE 


Equation  4-34  can  be  solved  at  any  temperature  for  any  desired  flux 
level  corresponding  to  the  given  operating  condition  of  the  detector, 
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P  =  P  (<t>,  T) 

The  detector  responsivity  is  given  by  4-2,  and  4-15 


(4-35) 


dP 

R  ~( - )<J>=>3 

34>  Y4> 


(4-36) 


where  <$>g  is  the  background  level  flux 


The  detector's  detectivity  is  then  proportional  to 


D* 


R  1/2 
(— ) 

P 


(4-37) 


The  results  which  will  be  displayed  below  have  been  obtained  using 
the  numerical  input  summarized  in  Table  4-1. 


TABLE  4-1.  INPUT  PARAMETERS  FOR  THE  CALCULATION  OF  THE  DETECTOR  RESPONSE 
(FLUX  4>B  =  0  or  4>B  =  1016  PR0T0NS/cm2-sec ) 


of 

Donor  Level 
the  Deep  Donor 

Acceptor  Level 
of  the  Deep  Donor 

Shallow  Acceptor 

Cnergy(eV) 

Eq=0 . 200 

EDa=0.75 

and  allowed  to  vary 

Ea=0 . 044 

Degeneracy 

y0 

90=(  M/2 

T+ 

9DA=(-~— )=2 

yO 

9A=( - )=4 

ta 

rhoton  Absorption 
Cross  Section 
(cm2) 

aD=5.4X10-16 

oDA=5.4x10~16 

aA=l  .8xl0-15 

Recombination 

Coefficient 

(cm3/sec ) 

BD=(6xlO~6/Ti -3) 

BDA=(6xlO-4/Tl -5) 

BA=(6xl0-4/Tl -5) 

Concentration 

N0=5xl016 

NA=5xl012 

( cm3) 
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The  system  will  be  investigated  for  a  variety  of  situations. 
First,  the  system  will  be  explored  in  thermal  equilibrium.  This  means 
the  light  will  be  off,  and  we  shall  study  both  high  and  low  counterdopant 
concentrations.  In  order  to  explore  the  consequences  of  the  negative  -U 
properties,  the  E  energy  will  be  allowed  to  vary  from  E  >  E 

U  M  U  M  U 

to  E  <  Eg.  The  same  properties  will  then  be  investigated  in  the 

presence  of  the  background  flux. 

Figure  4-3  provides  charge  center  concentrations  in  the  absence  of 

illumination  and  very  low  counterdopant  boron  concentration  (<j>  =  0, 

B 

Ng  =  5E12).  The  positive  -U  is  the  result.  E^  =  0.3eV  and 

E  =  0.225eV  show  that  the  deep  donor,  E^  =  0.2  eV,  exists  primarily 

in  its  neutral  charge  state.  Boron  is  totally  compensated  in  the  200K 
temperature  range.  Initially,  N+  =  Ng,  but  as  the  temperature 

rises  electrons  are  promoted  thermally  to  the  ( — /0 )  acceptor  level  and  a 
hole  is  created  in  the  valence  band.  This  transition  is  more  likely  than 
the  transformation  (0/+) 


D+  +  D°  +  hy  .  (4-38) 

This  is  true  because  there  are  many  more  than  N* 
(5  E  16  versus  5  E  12).  Therefore,  holes  are  created  predominately  via 
the  ( -/0 )  process. 


D°  -»  D~  +  hy  •  (4-39) 

Given  the  large  density  of  neutral  donors  still  surviving,  the  next 
process  is  the  capture  of  the  hole. 

D°  +  hv  -►  D+  .  (4-40) 

This  explains  why  we  eventually  see  for  a  positive  -U  in  Figure  4-3. 

Nfi  =  No  .  (4-41) 

These  two  charge  states  increase  in  tandem  at  the  expense  of  the  neutral 
donor  concentration.  At  the  same  time,  the  Fermi  level  locks  halfway 
between  the  (0/+)  and  (-/0)  energy  level.  (At  the  lowest  temperatures, 
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the  Fermi  level  is  locked  at  the  0/+  level.)  At  much  higher  temperatures 
thah  shown  in  Figure  4-3,  the  Fermi  level  is  unpinned  again,  N_ 

4-  U 

decreases,  increases,  and  Nq  decreases,  while  the  Fermi 

level  rises. 

The  behavior  of  the  negative  -U  center  system  in  Figure  4-3  is 

exemplified  by  the  0.15eV  (-/0)  and  0.20eV  (0/+)  system.  As  discussed 
earlier,  the  two  stable  charge  configuration  are  and  N*  with 

two  and  zero  electrons,  respectively.  The  introduction  of  neutral 

impurities  into  the  crystal,  at  low  temperatures,  produces  a 
thermodynamically  unstable  situation.  The  system  can  lower  its  energy  if 
half  the  impurities  lose  an  electron  and  the  other  half  pick  one  up,  i.e.; 

D°  -*  1/2  D+  +  1/2  D"  .  (4-42) 

(A  small  fraction  of  the  donors  give  up  one  electron  to  compensate  the 
residual  boron.)  Thus,  in  essence,  the  system  is  self-compensating. 


Nfj  =  N5  =  1/2  Nq  .  (4-43) 

This  results  in  the  Fermi  level  being  pinned  halfway  between  the  two 
energy  levels 


v  =  (Ed+Eda)/2  •  (4-44) 

As  the  temperature  increases,  electrons  from  the  valence  band  are 
promoted  thermally  to  the  first  available  empty  electronic  state  (0/+) 

D+  -»  D°  +  hy  ,  (4-45) 

since  putting  another  electron  on  D  is  ruled  out  (i.e.,  there  is  no 
(=/-)  state  in  the  gap) . 

Some  of  the  created  by  the  thermal  process  in  Equation  4-45  will 

+  — 

spontaneously  decay  as  in  Equation  4-42  into  D  and  D  .  The  thermal 
process,  Equation  4-45,  will  lead  to  an  increase  in  the  density  of 
neutral  donors  at  the  expense  of  the  concentration. 
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Another  process  that  leads  to  an  increase  of  the  D°  population  is 
the  hole  capture  by  0  . 


D“  +  hy  -»  D° 


(4-46) 


This  process  is  not  followed  by  rapid  hole 
at  low  temperature.  The  processes  given 
produce  the  increased  population  of 
Figure  4-3. 


emission  to  the  valence  band 
in  Equations  4-45  and  4-46 
centers  which  is  seen  in 


For  temperature  above  200K,  which  are  not  shown  in  Figure  4-3,  the 
o 

hole  captured  by  0  in  Equation  4-46  will  be  reemitted  to  the  valence 
band , 


D°  •*  D"  +  hy  .  (4-47) 

The  D  population  will  grow,  and  the  population  will  decrease.  At 
the  same  time,  there  will  be  a  rapid  decrease  of  the  D+  population 
since  the  thermal  ionization  of  D  leads  to 


D+  -*  0°  +  hy  -*•  D~  +  2hy  (4-48) 

In  the  second  step  of  this  process,  decays  via  the  process  of 

Equation  4-47  which  takes  less  energy  than  the  initial  ionization  of 

D  .  The  D  becomes  the  more  thermodynamically  stable  charge  state  at 

higher  temperature  because  the  Fermi  level  becomes  unpinned  again  and 

rises  above  y  given  by  Equation  4-44.  The  average  number  of  electrons 

per  donor  rises  from  "one"  at  low  temperatures  [N  =  1/2N  , 
0  u  u 

N5 =  0] 

navg  =  (N°  +2N0)/No  =  (2)(1/2)(ND/ND)  =  1  ,  (4-49) 

to  two  at  higher  temperatures  [Ng  =  Nq,  n£)  =  0] 

na vg  =  (2)(ND/N0)  =  2  .  (4-50) 
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4.8  THERMAL  EQUILIBRIUM  HIGH  BORON  CASE 


The  system  in  Figure  4-3  for  negative  -U  was  self-compensating  for 
lack  of  a  substantial  concentration  of  shallow  acceptors.  We  desire  to 
provide  a  comparison  in  Figure  4-4  with  the  system  in  Figure  4-3.  We  now 
make  the  density  of  shallow  acceptors  (Boron)  equal  to  that  of  the  donors 

NB  =  Nd  .  (4-51) 
As  the  donors  enter  the  silicon  matrix,  they  first  compensate  the  shallow 
acceptors . 

D°+A°  -*  D++A~  .  (4-52) 

The  equal  compensation  noted  in  Figure  4-3  means  that  there  are  no  more 
electrons  available  to  produce  the  auto  compensation  (Equation  4.1);  so 
at  low  temperatures,  we  have 


Nfi  =  N§  .  (4-53) 

The  Fermi  level  is  now  locked  half  way  between  the  boron  level  and  the 
(0/+)  acceptor  level  of  the  donor. 

This  same  fact  is  true  at  low  temperatures  for  the  positive  -U 
situation  in  Figure  4-4.  The  deep  donor  sim,  y  gives  up  its  electron  to 
the  shallower  boron  acceptor  and  the  Fermi  level  locks  halfway  in 
between.  As  the  temperature  rises,  the  conventional  process  occurs, 

D+  -♦  0°  +  hy  (4-54) 

D+  thermally  ionizes;  this  leads  to  an  increase  in  the  D® 
population.  As  the  temperature  continues  to  rise  thermal  ionization  of 

occurs, 

0°  -*■  D"  +  hy  (4-55) 

The  final  result  is  an  increased  population  of  the  D~. 
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Figure  4-4  Charge  Center  Concentrations  Ng  = 
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The  average  number  of  electrons  per  donor  increases  from  zero  at 
low  temperatures ,  as  the  Fermi  level  rises,  to  a  plateau  of  one  electron 
per  donor  as  the  Fermi  level  passes  the  (0/+)  level.  The  average  number 
of  electrons  reaches  a  value  of  two  per  donor  when  the  Fermi  level  rises 
above  the  (-/0)  level . 

The  temperature  behavior  of  the  negative  -U  system  in  Figure  4-4 
differs  substantially  from  that  of  the  positive  -U  system.  The  first 
step  is  to  thermally  ionize  the  (0/+)  level 

D+  -*  +  hy  (4-56) 

This  is  similar  to  the  case  for  the  positive  -U  system  in  Figure  4-4 

(Equation  4-54).  The  difference  is  in  the  subsequent  fate  of  some  of  the 

0 

neutral  donors.  The  next  transition  from  D  to  D  is  less  likely  to 

occur  in  the  positive  -U  system  than  in  the  negative  -U  system  at  the 

+  .  0 

temperature  where  0  is  easily  ionized  to  become  D  .  As  a  matter  of 
fact,  once  the  transition  in  Equation  4-56  is  accomplished  in  the 
negative  -U  system,  the  next  transition  required  even  less  energy 

D°  -»  D~  +  hv  (4-57) 

This  explains  why  the  D  population  grows  faster  and  the  D° 
population  slower  in  the  negative  -U  system.  The  average  number  of 
electrons  per  donor  atom  increases  rapidly  from  0  to  2,  without  a  plateau 
at  one  electron  per  donor,  as  the  Fermi  level  rises  with  temperature. 

4-9  BACKGROUND  ILLUMINATED  CASE 

An  understanding  of  carrier  dynamics  is  necessary  to  explain  the 
detector  performance  described  in  the  next  section.  In  a  typical  high 
background  application,  the  background  flux  is 

2 

4>  =  1  E  16  photons/cm  -sec 

D 
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The  carriers  dynamics  is  explored  using  the  same  equations  as  before,  but 
with  optical  generation  included.  Figures  4-5  and  4-6  detail  the 
behavior  of  the  low  and  high  counterdoped  systems,  respectively. 

4.9.1  Low  Boron,  with  Illumination 

The  presence  of  the  substantial  photon  flux  produces  a  large 
redi stribution  of  charge.  The  positive  -U  system  behavior  may  be 
analyzed  as  follows.  Essentially  all  donors  are  neutral  for  the  low 
temperature  -  no  light  case  in  Figure  4-3. 

n8  =  ND  . 

When  the  light  is  turned  on,  the  following  transition  occurs 

D°  +  hv  -*  O'  +  hv  (4-58) 

Simple  capture  of  the  carrier  produces 

D°  +  hv  -*  D+  ,  (4-59) 

D~  +  hv  -*  0°  .  (4-60) 

The  capture  process  by  the  neutral  donor  (Equation  4-59)  is  made  likely 
owing  to  the  large  concentration  of  neutral  donors;  this  is  in  spite  of 
the  low  recombination  cross-section,  Table  4-1.  The  system  is  self 
compensating  in  the  presence  of  light 

Nfi  =  Ng  (4-61) 

As  the  temperature  increases,  the  behavior  of  the  positive  -U 
system  depends  on  the  relative  positions  of  the  two  energy  levels.  It  is 
easier  to  reverse  the  process  in  Equations  4-58,  4-59,  and  4-60  (produce 
neutral  donors)  for  the  deeper  center  E  =  0.3, 

U  H 

(Thermal  Ionization)  D+  -*  D®  +  hy  ,  (4-62) 

(Capture)  D"  +  hy  -*  f)0  ,  (4-63) 
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(4-64) 


than  to  produce  further  thermal  ionization  of  0 

D°  -»  D~  +  hy 

This  leads  to  the  increase  of  N?  and  the  decrease  of  N+  and 

u  u  0 

N  At  still  higher  temperatures,  the  thermal  ionization  of  D 
becomes  likely  (Equation  4-64)  as  is  the  capture  of  a  hole. 

D®  +  hy  -*  D+  .  (4-65) 

So  that  Np  and  NQ  increase,  and  Np  decreases. 

0 

When  Eqa  =  0.225,  it  is  easy  to  thermally  ionize  D 

(Equation  4-64),  since  E^  ~  Ep  in  this  parametric  model.  This  is 

followed  by  the  capture  process  (Equation  4-65).  Therefore,  N+  and 

-  0  u 
Np  increase  at  the  expense  of  Np. 

The  explanation  for  the  behavior  of  the  negative  -U  system  is  as 
follows.  We  have  at  low  temperatures  in  thermal  equilibrium 

Nfi  =  Nq  =  1/2N0  , 

optical  absorption  leads  to 

D+  +  hv  ->  D° 

while  there  are  no  allowed  transitions  of  the  form 

D"  +  hv  ■*  D=  .  (4-68) 

What  happens  though  is  that  some  of  the  newly  created  neutral  donors 
decay  spontaneously  via 

D°  -»  D“  +  hv  ,  (4-69) 


(4-66) 

(4-67) 
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and  the  system  is  self-compensating  again.  The  overall  effect  of 
background  light  at  low  temperatures  is  the  production  of  an  otherwise 
metastable  population  at  the  expense  of  N*  and  NQ. 


As  the  temperature  increases,  the  population  is  thermally 


ionized , 


and  a  hole  is  captured 


D°  -+  0"  +  hy 


D°  +  hy  -♦  D+ 


(4-70) 


(4-71) 


At  still  higher  temperatures  the  (0/+)  level  comes  into  the  picture  and 


the  reverse  of  Equations  4-70  and  4-71  occurs, 

D+  -»  +  hy  (thermal  ionization  of  0/+) 


(4-72) 


D°  -»  D  +  hy  (thermal  ionization  of  ~/0) 


(4-73) 


This  replenishes  the  D°  population  at  the  expense  of  the  D+.  In 
Figure  4-5,  the  self-compensation  is  the  necessary  consequence  of  the 
very  low  boron  concentration  and  low  hole  concentration  in  the 
temperature  regime  of  interest. 

4.9.2.  High  Boron,  with  Illumination 

For  high  counterdopant  concentration,  the  thermal 
equilibrium  situation  is  depicted  in  Figure  4-4.  Essentially  all  the 
donors  are  used  up  to  compensate  boron 


Nj  =  ND  . 


(4-74) 


As  the  light  is  turned  on,  there  will  be  the  photoionization  process 

D+  +  hv  -*  D°  +  hy  ,  (4-75) 
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followed  by 


+  bv  -»  D  +  hy  .  (4-76) 

This  produces  a  substantial  Np,  NQ  density.  Most  holes  thus 

produced  are  captured  by  Ng  rendering  most  boron  neutral. 

Next,  we  discuss  the  positive  -U  system  as  the  temperature 

increases,  neutral  boron  begins  to  emit  holes  to  the  valence  band.  Those 

holes  that  are  recaptured  by  boron  are  even  faster  reemitted  to  the 

valence  band.  Therefore,  capture  at  the  boron  atoms  becomes  unlikely  as 

0 

the  temperature  increases.  The  capture  of  these  holes  by  0  and  0 

is  more  likely. 


D°  +  hy  -»  D+ 


(4-77) 


D"  +  hv  -»  D° 


(4-78) 


which  will  lead  to  a  decrease  of  the  D  and  D+  concentrations.  At 
still  higher  temperatures ,  the  thermal  ionization  of  the  (0/+)  level 
produces  a  hole. 


DJ  -»  D°  +  hy  (4-79) 

+  .  .0 

This  results  in  a  decrease  of  Np  and  an  increase  in  Np. 

The  negative  -U  systems  in  Figure  4-6  behave  similarly  for 
EpA  =  O.lSOeV  and  EpA  =  0.075eV.  The  original  populations  of  the 

charge  states  are  the  same  for  the  positive  and  negative  -U  systems  at 
low  temperatures  and  in  the  absence  of  light  (Figure  4-4).  The 
photoionization  of  the  D+  (Equation  4-75  and  4-75)  end  the  hole  capture 

by  N  shift  the  compensation  from  boron  when  th_  background  flux  is 

°  Ci¬ 

on.  Large  concentrations  of  D  ;  D  and  neutral  boron  are  created. 

As  before,  with  increasing  temperature  it  is  the  turn  of  boron  to  emit 

holes  to  the  valence  band.  These  holes  are  captured  by  0^  and  D  , 

(Equations  4-77  and  4-78)  as  in  the  case  of  the  positive  -U  system  in 
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Figure  4-6.  The  EQA  =  0.075eV  sy  tem  shows  that  with  a  further 
increase  of  temperature  the  (-/0)  level  participates  in  the  thermal 
ionization 

D°  -»  D_  +  hy  .  (4-80) 

This  produces  an  increase  in  Nq  and  a  decrease  in  N^. 

A  further  increase  in  temperature  involves  the  (0/+)  level  in  the 
conversion  D+  to  D  which  thus  increases  at  the  expense  of 

D+  ->  D°  +  hy  (4-81) 

Of  course,  some  of  the  D^1  will  be  converted  to  D  (Equation  4-8C) 

since  E„.  <  E„.  This  is  the  main  difference  between  the  positive 
DA  0 

and  negative  -U  systems  in  Figure  4-6. 

4.10  RESPONSI VITY 


Responsivity  (Equation  4-36)  is  a  measure  of  the  detector's  ability 
to  detect  incident  signal  radiation  while  the  detector  is  at  temperature 
"T"  and  is  already  exposed  to  the  background  radiation.  The  partial 
derivative  in  Equation  4-36  measures  'he  additional  photohole  density 
Ap  when  additional  flux  "A4>"  falls  on  the  detector. 

Figure  4-7  displays  the  calculated  responsi vities  for  several 
systems  of  interest.  The  behavior  of  responsivity  is  easily  understood 
in  the  high  temperature  region  T  >  100  K.  We  recall  that  the  principal 
photoconducti ve  transition  is  given  by 

D+  +  hv  -»  D°  +  hv  '  (4-82) 

This  results  in  high  quantum  efficiencies  and  responsi vities  being 
promoted  by  large  D+  concentrations.  Figure  4-6,  for  high  dopant 
concentrations,  shows  that  increases  with  temperature  as  the 
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Figure  w  Calculated  Responsivity 


r  i 


compensation  is  shifted  back  to  boron.  This  makes  the  material  more 

sensitive  to  the  incoming  radiations,  which  means  it  is  more  responsive. 
The  D+  are  thermally  ionized  above  100K;  this  produces  and 

results  in  a  drop  in  responsi vity . 

The  negative  -U  system  =  0.075eV  has  about  five  times  as  high 
a  responsivity  as  the  sy  .terns  with  E  >  0.15.  This,  the  most 

U  n 

negative  -U  system  considered  in  this  work,  has  a  more  efficient  transfer 

0 

of  holes  from  boron  to  D  (Equation  4-77), 

D°  +  hy  ->  D+  (capture)  (4-77) 

It  also  is  most  efficient  in  thermal  emission  of  holes  by  because 
Eda  i s  small . 

-»  D“  +  hy  (emission)  (4-80) 

Capture  of  the  hole  is  the  next  step.  This  gives  rise  to  an  enhanced 
D+  population  which  is  responsible  for  the  large  responsivity. 

The  counterdoping  concept,  with  the  large  boron  concentration, 
works  because  donors  compensate  boron.  Thus  the  desired  deep  donor  level 
is  exposed.  This  works  for  both  the  positive  and  negative  -U  systems. 
What  is  not  appreciated  is  that  the  counterdoping  scheme  works  without 
any  boron  if  the  center  has  the  negative  -U  property.  As  discussed  in 
connection  with  Figures  4-3  and  4-5,  for  low  boron  concentrations,  the 
system  is  self-compensating.  The  responsivity  for  the  En.  =  0.075eV 

3  r* 

( Ng  =  5E12  cm  )  system  in  Figure  4-7  is  within  an  order  of  magnitude 
of  the  high  boron  case.  This  disadvantage  may  be  offset  by  what  is 
expected  to  be  the  larger  hole  mobility  in  the  boron  free  negative  -U 
system. 

4.11  DETECTIVITY 


The  thermal  generation-recombination  (G-R)  noise  limits  the 
detectivity  of  photoconducti ve  infrared  detectors.  Detectivity  is 
therefore  degraded  when  there  is  a  large  scale  thermal  emission  of  holes 
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to  the  valence  band.  The  detector  is  background-limited,  at  low 
tempci  atures,  when  the  only  aperture  noise  mechanism  is  the  baui<.y  round 
photon  induced  generation-recombination  noise  (Figure  4-8).  The 

detectivity  is  then  temperature  independent.  The  first  temperature 
roll-off  in  the  detectivity  occurs  when  boron  is  thermally  ionized  for 
the  high  boron  case.  The  large  increase  in  the  hole  population  gives 
rise  to  the  G-R  noise.  These  holes  are  subsequently  recaptured  by  the 
donors  and  the  detectivity  flattens  out.  At  higher  temperatures,  the 
emission  of  holes  from  the  (0/+)  donor  level  produces  the  final 
temperature  roll-off  for  T  >  TOOK. 


The  low  boron  case,  E^  =  0.075eV, 

first  roll-off  at  higher  temperature 

E  =  0.044eV.  Here,  the  first  thermal 
A 

emission  by  the  neutral  donors  with  the  E 
systems  compared  in  Figure  4-8  show  compara 


negative  -U  case,  shows 

than  the  high  boron 

roll-off  is  due  to  the 

activation  energy.  All 
DA 

ble  detectivities. 


the 

case 

hole 

the 
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5.0  CONCLUSIONS 


The  counterdoping  scheme  has  been  modeled  for  systems  exhibiting 
multivalent  statistics.  Interesting  differences  were  observed  for 
systems  with  negative  and  positive  -U  properties.  This  may  well  be  the 
first  such  investigation  of  the  infrared  detector  potential  of  the 
negative  -U  center  and  of  the  multivalent  statistics  in  the  counterdoped 
scheme . 


We  find  that  positive  and  negative  -U  systems  show  similar 
detectivities  and  T^  temperatures .  Most  interesting,  for  detector 
applications,  is  that  the  negative  -U  systems  are  self-compensating,  or 
sel f-counterdoped ,  because  of  the  inverted  level  ordering.  They  appear 
to  have  comparable  responsi vities  to  that  of  the  boron  counterdoped 
detectors.  More  quantitative  predictions  await  better  experimental  input 
to  the  theory. 

The  model  outlined  in  the  theoretical  section  has  identified  the 
relevant  parameters  which  control  the  operating  temperature  and  figure  of 
merit  for  the  counterdoped  detector.  The  rate  equation  can  be  used  to 
predict  the  response  of  any  counterdoped  system.  In  particular,  it  can 
be  used  to  optimize  doping  concentrations  for  the  highest  BLIP  operating 
temperature . 

We  have  successfully  created  the  positive  divacancy  and  observed  to 

presence  using  DLTS ,  K'll  and  infrared  analysis.  The  presence  of  boron 

and  carbon  in  the  silicon  sample  lead  to  creation  of  other  complexes 

which  were  also  observed.  These  are  the  interstitial  carbon 

substitutional  carbon  trap  at  E  +  0.44eV  which  was  observed  in  the 

v 

DLTS  spectra  [Section  3.1  and  3.2]  as  well  as  the  boron  related  trap 
which  closely  overlaps  the  positive  divacancy. 

The  final  statement  is  that  the  evidence  produced  in  the  infrared 
spectral  measurements  shows  that  illumination  does  not  upset  the  creation 
of  positive  divacancies.  The  counterdoping  concept  is  valid  in  that 
boron  compensation  serves  to  create  positive  divacancies  and  prevents 
formation  of  neutral  and  negative  divacancies. 
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ABSTRACT 

A  parametric  study  of  a  counterdoped  extrinsic  silicon 
detector  system  is  presented.  The  system  consists  of  a  shallow 
group  III  acceptor  and  a  deep  donor  with  a  level  near  the  valence 
band.  It  is  known  that  compensated  donors  become  neutral  upon 
photoabsorption  and  present  a  very  small  capture  cross  section 
for  the  carriers.  This  leads  to  a  higher  BLIP  operating  temper¬ 
ature  for  the  counterdoped  system  than  for  the  conventional 
detector  operating  at  the  same  wavelength.  A  detailed  micro¬ 
scopic  model  for  the  counterdoped  system  is  developed  to  study 
the  temperature  dependence  of  responsivity  and  detectivity. 
Detectivity  is  studied  as  a  function  of  the  hole-neutral  donor 
recombination  coefficient  and  as  a  function  of  the  acceptor-to- 
do.nor  density  ratio.  The  best  detectivities  and,  simultaneously, 
responsivities  are  obtained  for  acceptor  density  equal  to  that 
of  the  donors.  Using  boron,  aluminum,  gallium,  or  indium  as 
the  group  III  acceptor  in  the  parametric  study  shows  that  the 
highest  responsivities  and  detectivities  are  obtained  for  the 
shallowest  acceptor,  that  is,  boron.  These  results  can  be  under¬ 
stood  on  the  basis  of  the  picture  involving  optical  and  thermal 
carrier  generation,  carrier  recombination,  and  the  dynamic  inter¬ 
action  of  the  donor  and  acceptor  levels  mediated  by  the  valence 
band . 


1.0  INTRODUCTION 

Extrinsic  silicon  infrared  detectors  in  the  3-5  urn  range  are  limited  to  oper¬ 
ating  temperatures  lower  than  60K.  High  BLIP  (Background  Limited  Infrared  Photo¬ 
conductor)  temperature  silicon  infrared  detectors  are  required  to  reduce  the  size 
and  weight  of  the  cooling  systems  associated  with  the  detector  arrays.  The  limita¬ 
tion  or.  the  detector  operating  temperature  can  be  traced  to  the  large  recombina¬ 
tion  coefficients  of  carriers  in  conventionally  doped  extrinsic  silicon. 

A  scheme  for  increasing  the  BLIP  operating  temperature  of  an  extrinsic 
detector  has  been  proposed  by  Elliott,  Migliorato,  and  Vere'*'.  Specifically,  they 
proposed  to  dope  a  semiconductor  with  either  a  deep  donor  or  a  deep  acceptor  and 
then  counterdope  the  material  with  impurity  of  opposite  kind.  The  principal  IR- 
active  centers  are  the  deep  impurities  which  have  lost  charge  to  the  compensating 

*  This  work  is  supported  by  the  Materials  Laboratory,  Air  Force  Wright  Aeronautical 
Laboratories,  Contract  No.  F3 36 1 5- 8 5-C-5062 . 

Elliott,  C.  T.,  Migliorato,  P.,  and  Vere ,  A.  W. ,  "Counterdoped  extrinsic  silicon 
infrared  detectors,"  Infrared  Physics  18,  63  (1978). 
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counterdopant.  Therefore,  the  deep  acceptors  become  negatively  charged  and  deep 
donors  become  positively  charged. 

In  the  conventionally  doped  extrinsic  detectors  IR  excitation  produces  photo¬ 
carriers  of  one  sign  and  an  immobile  center  of  opposite  sign.  The  attractive 
interaction  between  the  carriers  and  the  photoionized  centers  leads  to  short 
carrier  lifetimes  which  ultimately  results  in  low  detector  operating  temperatures. 
On  the  other  hand,  in  the  counterdoped  detector  the  photoionized  center  becomes 
neutralized  in  the  photoabsorption  and  can  be  either  neutral  or  even  repulsive  to 
the  created  photocarriers.  As  the  result  the  recombination  coefficients  are 
decreased  and  carrier  lifetimes  are  increased,  which  leads  to  higher  BLIP  oper¬ 
ating  temperatures. 

The  original  proposal  of  Elliott,  Migliorato,  and  Vere  contains  a  theoretical 
analysis  of  the  counterdoped  system.  In  particular,  these  authors  provided  expres¬ 
sions  for  spectral  responsivities  and  detectivities  as  well  as  for  the  BLIP  tem¬ 
perature.  The  present  work  will  provide  a  more  detailed  microscopic  model  of  the 
counterdoped  system  and  will  explore  the  detector  performance  parameters  over  the 
whole  temperature  range.  By  examining  the  carrier  dynamics  as  a  function  of 
doping  conditions  and  recombination  coef f icients ,  specific  recommendations  can  be 
made  for  optimum  choice  of  detector  parameters.  It  is  hoped  that  the  present 
model  will  add  to  the  understanding  of  the  physical  basis  for  the  counterdoped 
detector . 


Section  2  describes  the  model  of  the  counterdoped  system,  while  Section  3  is 
devoted  to  the  presentation  of  results.  Section  4  contains  the  principal  conclu¬ 
sions  of  this  work. 


2.0  MODEL  OF  THE  COUNTERDOPED  SYSTEM 
2.1  DESCRIPTION  OF  THE  DETECTOR 

Figure  1  displays  the  schematic  band  diagram  for  the  counterdoped  semiconduc- 
toi .  The  principal  infrared  active  dopant  is  the  deep  donor  with  energy  level  in 
the  gap  at  energy  E^  from  the  valence  band.  The  donor  concentration  is  N^.  For 
the  sake  of  simplicity,  the  center  responsible  for  the  donor  formation  can  only  be 
either  neutral  or  singly  positively  charged.  Extension  of  the  present  formalism 
to  multivalent  flaws  or  impurities  should  be  straightforward. 


The  system  is  also  doped  with  acceptors  with  density  at  energy  Eft  above 
the  valence  band  corresponding  to  the  0/-  charge  transition  of  the  center.  The 
simplest  and  in  practice  best  realization  of  the  acceptors  are  the  simple  hydro- 
genic  column  III  impurities:  boron,  aluminum,  gallium,  or  indium.  At  low  temper¬ 
atures  donors  lose  one  electron  each  to  the  acceptors  below  so  that 


or 


for  nd  I  na 


N .  for  N„  >  N„ 
A  D  —  A 


(1) 


(2) 


Figure  1.  Schematic  energy  level  system  for  the  counterdoped  system  containing 
a  shallow  acceptor  and  a  deep  donor  level.  EG ,  Ev ,  En»  and  EA  stand 
for  the  conduction  band,  valence  band,  the  donor  level,  and  acceptor 
level  energy  positions,  respectively. 

Photoconductivity  results  from  photoabsorption  of  a  photon  with  energy 
hv>_ED#  by  a  valence  electron  which  is  then  promoted  to  the  donor  lev  1 

n+  +  hv  -  D°  +  hy  .  (3) 

As  the  result  the  donor  is  neutralized  and  a  hole  in  the  valence  band,  hy,  is 
created.  The  photocurrent  is  carried  by  the  valence  band  holes.  Equivalently, 
of  course,  one  may  think  of  Eq.  (3)  as  a  promotion  of  the  hole  on  the  donor  level, 
D+ ,  to  the  valence  band. 

The  hole  lifetime  in  the  valence  band,  as  the  hole  is  swept  out  by  the  elec¬ 
tric  field  in  the  detector,  ought  to  be  long  since  the  recombination  cross  section 
for  the  hole  capture  by  neutral  acceptors  should  be  very  small.  The  recombination 
process  would  proceed  via  the  reaction 

D°  +  hy  -  D+  +  (energy)  (4) 

which  has  to  be  contrasted  with  the  recombination  process  in  the  standard  p-type 
photoconductive  detector  where 

A  +  hy  -  A°  +  energy  .  (5) 

In  Eq .  (5)  the  recombination  cross  section  for  the  hole  capture  by  the  attractive 

center  is  estimated  to  be  several  orders  of  magnitude  larger  than  for  the  process, 
Eq .  (4),  in  the  counterdoped  detectors.  For  multivalent  flaws  it  is  possible  to 

have  the  capturing  center  to  be  of  the  same  electrical  sign  as  the  photocarrier, 
which  would  further  reduce  the  recombination  cross  section. 

In  order  to  calculate  the  two  f igures-of-merit  for  the  photodetector,  the 

responsivity  and  detectivity,  the  rate  equations  for  the  optical  and  thermal 

2 

carrier  generation  must  be  solved  .  The  rate  at  which  neutral  acceptors  are  pro¬ 
duced  is  given  by 

^Szmulowicz,  Frank,  Hemenger,  Patrick,  and  McGuigan,  Shaun,  "Low  background  per¬ 
formance  modeling  of  gallium  doped  silicon  infrared  detectors,"  Proc .  Materials 
IRIS  Conference,  Colorado  Springs,  CO  (1985). 


(6) 


where 


3Na  „  _  B, 

TT  -  pBANA  -  W*  -  NA  rk  NV  exp('VkBT) 


c , 

e\ 


c 


the  degeneracy  of  the  acceptor  level,  g^=4 

the  recombination  coefficients  for  holes  with  the  ionized  acceptor,  Eq .  5 
valence  band  hole  concentration 

N.  =  total  acceptor,  neutral  acceptor,  and  ionized  acceptor  concentra¬ 
tions,  respectively 

wavelength  dependent  photoabsorption  cross  section  for  the  acceptor  level 
incident  photon  flux  consisting  of  the  background,  <>B,  and  signal,  ps, 
radiation 

convolution  of  the  photon  flux  with  the  photoabsorption  cross  section 
effective  valence  band  density  of  states 
acceptor  photoionization  energy  . 


The  convolution  <ct>  will  be  replaced  by  c  *,  where  c  is  the  peak  photoabsorption 
"  " "  p  P 

cross  section  for  the  given  impurity.  The  replacement  thus  defines  an  effective 

flux  t. 


The  first  term  in  Eq .  (6)  represents  the  hole  capture  rate  by  ionized  accept- 

-4  1.5 

ors.  In  the  present  parametric  study,  3.  will  be  taken  to  be  E  =6x10  /T 
(cm.  /sec) .  The  temperature  dependence  of  B  is  consistent  with  the  results  of 
Geim,  Pensl,  and  Schultz3  on  indium  and  gallium  in  silicon  and  with  the  theore¬ 
tical  model  of  Lax4.  The  numerical  prefactor  is  also  typical  of  the  results  of 
Geim  et  al.  The  second  term  in  Eq .  (6)  gives  the  photoabsorption  rates  of  holes 

from  neutral  acceptors  and  the  third  term  gives  the  corresponding  thermal  hole 
generation  rate.  The  typical  peak  photoabsorption  cross  section  for  simple 
acceptors  in  silicon  is  in  the  c  =10  33-10  3^cm3  range.  Given  that  in  steady  state 


3N,/3t  =  0,  and 


(7) 


gives 


where 


aANA 

N  = 

A  p+a. 


0  A  =  (Nv/9A)e>Cp(_EA/kBT:  +  <0A*VBA 


(8) 

(9) 


For  donors  the  rate  equation  is  given  by 


3N. 


aT~  =  bdpnd  -  gDBDNvexp(-EDABT,ND  -  nd<cd4> 


(10) 


where 

N’  ,N°,N*  =  total,  neutral,  and  ionized  donor  densities 
DUD 

BD  =  the  recombination  coefficient  for  holes  with  neutral  donors  (to  be 
varied  in  the  calculation) 


^Geirrw  K .  ,  Pensl ,  G.,  and  Schultz,  M.,  "Shallow  acceptor  population  and  free  hole 
concentration  in  Si:In  and  Si :Ga  with  IR-photoexcitation , "  Appl.  Phys.  A2_7, 

71  (1982). 

4Lax,  M. ,  "Cascade  capture  of  electrons  in  solids,"  Phys.  Rev.  119 ,  1502  (1960). 


g3  «  the  degeneracy  factor  for  the  donor  leve*,  *aken  to  be  9D=2  for  the 

sake  of  concreteness 

Ep  =  the  donor  level  position  with  respect  to  the  valenc=  band 

cQ  =  the  photoab.orption  cross  fiction  of  the  ionized  donor  level,  Eq .  (3). 

Please  note  the  different  placement  of  the  degeneracy  factors  in  Eqs.  (6)  and 
(10).  These  can  be  shown  to  be  necessitated  by  the  detailed  balance  conditions. 
The  first  term  in  Eg.  (10)  is  the  hole  capture  rate  by  neutral  donors,  the  second 
and  third  are  the  thermal  and  optical  excitation  rates  of  holes  to  the  valence 
band  from  the  ionized  donors,  Eq .  (3).  A.~ain,  in  the  steadv  state  3.\'*/3t=0,  and 


Eq .  ( 10 )  i  ields 


o  + 
ND  +  ND 


<  =  — 

D  aD+p 


an  =  N,.gr,exp(-En/k„T)  +  <cn*>/Br 


Using  the  charge  neutrality  condition 


P  =  N, 


together  with  Eqs.  (8)  and  (12),  result  in  the  cubic  equation  for  p, 

p3  +  p2[ND-aA+aD]  +  oaA[aD+ND-NA]  -  =  0  •  (15) 

Equation  (15)  has  to  be  solved  for  p  at  different  temperatures  and  background 
flux  levels.  In  special  circumstances  Eq.  (15)  reduces  to  an  easily  solveable 
quadratic  or  even  a  linear  equation.  In  genera,,  though,  one  must  solve  the  cubic 
equation . 

2.2  DETECTOR  PERFORMANCE  MODEL 

Spectral  responsivity  is  defined  as  the  signal  current  per  unit  incident  sig¬ 
nal  radiation  power.  The  signal  current  is  proportions1  to  th<_  density  of  signal 
generated  carriers  pg  which  is  the  difference  between  the  densitie.  of  holes  with 
and  without  the  signal  flux  on,  i.e. 

P5  =  PUB+4S.T)  -  p(i~,T)  .  (16) 

In  the  small  signal  approximation  we  can  take,  from  Eq .  (If) ,  the  signal  flux 
generated  carrier  density  to  be 

Ps  ( 3p/3 <t )  ^  $s  .  (17) 

B 

Since  the  signal  power  is  also  proportional  to  $  ,  the  responsivity  is  proportional 
to 

R  v  (3p/34)  (18) 


with  the  derivative  evaluated  at  $=fri. 

ts 

The  spectral  aetectivity  is  proportional  to  the  ratio  of  the  resp^"  ivity 
to  the  root-mean  square  noise  cur: ent 

V*  v  R/I..  .  as 


10  20  30  SO  70  100  Sco 

'^V?r^T VRZ  (X) 

Figure  4-7  Calculated  Responsivity 
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The  genera t ion-recombinat ion  noise  current  is  given  by 

XN  "  (4Ioe  Gif)l5  '  (20> 

where  IQ  is  the  background  photon  and  thermally  generated  current 

x0  '  P(eB.T)  . 

Li  is  the  band width  of  the  measuring  electronic  circuit,  and  G  is  the  gain  of  the 
detector.  The  gain  is  the  ratio  of  the  number  of  holes  per  second  flowing  througi 
the  detector  per  photon  absorbed  in  it,  i.e.. 


;/  3  i  i  .  ■*  ,  from  Eq 
*  3  s 


?r ,  Eqs.  (20-2 2)  give 


G  -v  Is/0s  . 

(17),  the  gain  becomes 

G  “V  (  op/ 5$)  ^ 

’B 

I  fip.'Jo) -  = 

V — s— ^  • 


E-;s.  (18)  and 

cf  merit  for 
to  display  a: 


(15).  The  ce: 


(21) 


(22) 


(23) 

es 


(23)  constitute  the  primary  means  for  evaluating  the  two  figu 
ihe  detector.  In  the  present  parametric  study  it  is  less  importer, 
solute  magnitudes  of  the  performance  figures  than  it  is  to  compare 
:u  of  detectors  as  various  detector  parameters  are  varied. 

density  p  in  Eq .  (23)  is  calculated  from  the  cubic  equation,  Eq . 

nvative  cf  p  with  respect  to  the  photon  flux  is  given  by 
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ndp 
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whicl  is  symmetric  in  A,0  indices. 

3.0  RESULTS 


figure  2  shows  the  calculated  detectivity  for  the  input  parameters  of  Table 
Not*..*  tout  the  donor  level  has  been  placed  arbitrarily  at  the  energy  of  gallium, 
acceptor  ir.  silicon,  with  the  donor  concentration  and  its  absorption  cross  sect  i  c 
also  typical  of  gallium.  Therefore,  Figure  2  can  also  be  used  to  compare  a 
counter d_-p**d  and  a  conventional  detector  operating  in  the  same  wavelength  ranee. 

I  r.  Fiiure  2  the  comp^nsat  i  o.n  is  exact  N.-N-,  and  the  recombination  coefficient  is 

M  b 

varied  ov-r  orders  of  m.a  :r.  i  t  ude  from  that  typical  of  attractive  centers  in 

c*.:r.  1 1  ■.  n.j  1 1  y  d  d  (detectors  t  the  small  recomb  i  na  1 1  or.  coefficients  expected  c 

n*:-j *  r 3  1  cr  repu  1  s i v-*  centers  in  ecus tt»rd oped  detectors.  Enina  the  input  ir.  Table 
1,  Figure  2  i  n  \ :  \  *  ■  ■  a  the*  possibility  of  rain  mi  T^.,p  from  about  25  K  to  5  5b  for 

th>*  detect  or  oj. ‘rating  m  the1  8-13  wncbw.  I  hasten  to  add  that  at  present 
r*»  is  r  ■  ■  s.;*.ible  d* '  r  d->nor  in  this  wavelength  ramie . 
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Figure  2.  The  detectivity  of  a  model  counterdoped  system  as  a  function  cf  the 
recombination  coefficient  and  temperature.  The  detectivity  is  dis¬ 
played  on  a  relative  scale. 

Table  1.  Input  parameters  for  the  results  of  the  calculation  displayed  in 
Figures  2,  3,  and  4.  The  background  level  is  taken  to  be 
vn  =  10^  photons/cm2-sec . 
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B (cm  /sec) 

Vari able 

6x10~4/T1 ' 5 

concentration  (cm  3) 

5xl016 
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energy  (meV) 

72 

45 

The  choice  of  exact  compensation  for  the  display  in  Figure  2  and  the  shap 
of  the  detectivity  curve  has  to  be  explained.  In  Figure  3  BD  is  fixed  at 
6xl0"10/T1 - 5  (cm3/sec)  and  at  5x10 1 6  cm” 3 ,  while  the  acceptor  compensating 
density  is  varied.  The  choice  of  the  NA/ND  ratio  will  have  to  be  deferred  unt 
corresponding  responsivity  curves.  Figure  4,  are  discussed. 

Let  us  begin  the  discussion  of  Figure  3  with  the  case.  At  the  lr»e 

of  temperatures,  without  the  background  flux,  all  the  donors  and  acceptors  are 
ion l zed ; 
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Figure  3.  The  detectivity  (on  a  relative  scale)  of  a  model  cour.terdoped  system 
as  a  function  of  acceptor  concentration  and  temperature  for  a  fixed 
donor  concentration  and  B  =6xl0~  10/T1  • 5  (cmVsec)  . 


As  the  background  flux  is  turned  on  optical  absorption  of  valence  electrons  can 
proceed  to  the  donor  level  only, 

D+  +  hv  -»  D°  +  fly  ,  (26) 

since  the  acceptor  levels  are  already  occupied  by  electrons.  The  valence  hole  is 
more  likely  to  recombine  at  an  ionized  acceptor  A  , 

A  +  h,^  =  A°  +  (energy)  ,  (27) 

than  at  the  neutral  donor  because  of  the  widely  different  capture  cross  sections. 
Therefore,  in  Figure  3  for  T> 20K  with  the  background  flux  on, 
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(28) 

(29) 


a .  =  c.f/B.  <■  «■  p  <  <  N,,  N„  =  a.N./p  , 

A  A  A  A  A  /A  A 

aD  =  CD*/  BD  ' J  P'  >:>  P  '  ND  =  NDP'/;ir 
and  the  charge  neutrality  equation  is  given  by 

p2  [ND/aD  +  1]  =  aANA  . 


The  solution  for  p  becomes 

so  that. 


(30) 

(31) 

(32) 

(33) 

(34' 


independent  of  temperature.  This  simply  means  that  the  detector  is  background 
limited  before  the  onset  of  significant  thermal  hole  generation. 


The  optical  neutralization  of  acceptors  and  donors,  Eqs.  (28-29),  and  the 
lack  of  thermally  generated  carriers  explain  the  flat  portion  of  the  V‘  vs  T  curve 
for  T<20K.  As  the  temperature  is  increased  holes  are  thermally  ionized  from  the 
boron  acceptors  to  the  valence  band, 

A°  -  A-  +  hy  ,  (35) 

which  contributes  to  the  thermal  g-r  noise.  Holes  in  the  valence  band  may  be 
recaptured  at  A  centers,  but  then  are  rapidly  thermally  reemitted  to  the  valence 
band.  Therefore,  a  more  likely  process  is  the  capture  at  the  previously  neutral 
donors 

D°  +  hy  -  D+  .  (36) 

This  process  of  hole  capture  by  donors  from  acceptors  mediated  by  the  valence  band 
lowers  the  detectivity  because  of  the  thermal  rolloff  owing  to  the  thermal 
ionization  of  boron.  At  the  same  time,  this  process  will  be  seen  to  improve  the 
responsiv ty  since  the  concentration  of  the  D  principal  IR-active  level  increases 

Analytically,  the  region  above  T>20X  is  characterized  by  the  following 
conditions , 

aA  "  P  “  NA  '  NA  =  NaVP  '  (37) 

Nd  >>  aD  >>  P  ,  N*  =  NDp/aD  ,  (38) 

as  before,  except  that  now  cA  possesses  a  significant  thermal  component.  The 
solution  for  p  is  again  given  by 

p  =  ( ) **  '  (39) 

and  V*  is  given  by 

^[VVd1  +  V(Va’^  •  (40) 

The  two  terms  in  Eq.  (39)  are  of  about  the  same  order  of  magnitude  for  20-T>-25K, 
except  that  increases  rapidly,  owing  to  the  thermal  depopulation  of  the  boron 
level,  while  is  relatively  constant.  This  then  provides  the  explanation  for 
the  thermal  rolloff  starting  at  T=20K  in  Figure  3. 


fhe  holes  thermally  emitted  by  the  boron  level  experience  different  fates 
ling  on  the  relative  concentrations  of  boron  with  respect  to  doners.  For 
N  the  holes  are  quickly  captured  by  the  donors  so  that 


an  :  the  thermal  rolloff  in  the  detectivity  is  stopped.  As  can  bo  seen  in  Figure 
curves  1,  2,  and  3  tend  back  toward  the  background  limited  behavior. 


For  N.  -  S'  the  transfer  of  holes  from  the  boron  level  to  the  donor  level  is 
not  complete  at  low  temperatures.  Therefore,  a  substantial  density  of  neutral 
boron  remains  which  contributes  to  the  thermal  noise.  In  particular,  curves  5 
and  6  in  Figure  3  with  exhibit  lower  detectivities  for  just  that  reason. 


ture  is  raised  there  is  also  a  thermal  depopu 
'he  thermal  rolloff  due  to  the  donors  occurs  a 
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ure  3  exhibit  comparable  detectivities,  but  t 
ary  by  orders  of  magnitude.  Several  responsi 
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carrier  lifetime.  Increasing  the  temperature  further  results  in  thermal  depopula¬ 
tion  of  the  donor  levels  so  that  the  quantum  efficiency  and  lifetime  decrease  as 
does  the  responsivity . 

The  principal  conclusions  that  emerge  from  this  part  of  the  parametric  model¬ 
ing  effort  are  as  follows.  First,  a  dramatic  improvement  in  the  Tn,T_  is  realized 
for  counterdoped  systems  over  otherwise  comparacle  conventionally  doped  systems. 
For  the  single-valent  donor  system  the  best  responsivities  and  simultaneously 
best  detectivities,  with  the  the  highest  are  achieved  for  equal  compensa¬ 

tion  of  the  donor  concentration. 

With  these  preliminary  considerations  out  of  the  way,  we  may  turn  to  a  more 
realizable  system8  This  system  consists  of  the  positive  divacancy  level  0/+ 

at  E  =200  neV  with,  the  peak  resoonse  at  4  pm.  The  vt  oositive  divacancies  can  be 

J  '  2  * 

produced  by  e-beam  irradiation  of  p-type  silicon  with  the  dose  adjusted  to  match 
the  density  of  acceptors.  Proper  annealing  procedure  must  be  followed  to  repair 
the  lattice  damage  while  maintaining  a  high  density  of  the  divacancies.  I  model 
this  system  using  the  parameters  from  Table  1,  except  as  noted  below.  The  peak 
photoabscrption  cross  section  for  the  +/0  divacancy  transition  is  estimated  to  be 
4x10  ^cr.2  by  Humphreys  et  al8  and  2.5x10  2‘6cm2  by  Brotherton  et  al8  which  is  con¬ 
sistent  with  the  value  in  Table  1  adopted  for  the  divacancy  in  the  present  work. 
For  modeling  purposes,  I  adopt  B„=6xl0  ^/T1'8  cm'Vsec,  although  the  temperature 
dependence  of  the  hole  recombination  coefficient  at  neutral  divacancy  is  not 
known.  The  numerical  constant  has  been  adjusted  to  correspond  to  the  measured 
recombination  cross  section  in  the  71-95K  range  of  4x10  16cm2  obtained  by 
Brotherton  et  al  via  DLTS.  The  donor  activation  energy  is  taken  to  be  ED=200  meV. 

There  is  a  choice  of  acceptors  for  the  p-type  dopant  of  the  starting  mater¬ 
ial:  boron,  aluminum,  gallium,  and  indium.  Figure  5  shows  the  detectivity  and 

responsivity  for  these  four  choices  of  the  starting  p-type  material.  The  respon¬ 
sivity  for  indium  doping  is  too  small  to  be  displayed  on  the  scale  of  the  drawing. 
From,  the  detectivity  curves  it  appears  that  indium  might  be  the  preferred  choice 
given  its  high  thermal  activation  energy  which  delays  its  thermal  rolloff  tj 
higher  temperatures  than  for  the  other  p-type  dopants.  The  low  value  of  the 
responsivity  using  indium  doped  silicon  mitigates  against  this  choice. 

The  problem  with  indium  can  be  explained  by  tie  following  argument.  As 
explained  earlier,  responsivity  increases  with  the  increase  in  the  D+  concentra¬ 
tion.  D+  are  created  by  the  hole  capture  from  the  valence  band  with  the  holes 
themselves  thermally  emitted  by  the  acceptors.  The  deeper  acceptors  start  emit¬ 
ting  holes  thermally  at  higher  temperatures  than  does  the  shallow  boron  acceptor. 
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Figure  5.  Detectivity  and  responsivity  (on  a  relative  scale)  of  a  counterdoped 
detector  with  the  divacancy  deep  donors  as  the  principal  IP.-level  and 
with  four  choices  for  the  acceptor  counterdopant:  boron,  aluminum , 
gallium,  and  indium. 

In  particular,  at  the  temperature  where  there  is  a  significant  transfer  of  holes 
from  indiun  to  the  divacancy,  the  divacancy  itself  begins  to  emit  holes  thermally 
back  to  the  valence  band.  On  the  other  hand,  boron  is  well  into  the  exhaustion 
region  with  all  the  holes  captured  by  the  divacancies  before  the  divacancy  caused 
thermal  rolloff  begins  at  about  90K.  This  is  consistent  with  the  TgT  of  100K, 
reported  by  Maher  et  al,  for  a  similar  background  flux. 

Figure  5  indicates  that  the  best  choice  for  the  p-type  starting  material  is 
that  containing  the  shallowest  acceptor.  It  is  important  that  the  knee  of  the  F* 
curve,  which  determines  T„it_,  is  at  the  same  temperature  as  the  peak  in  the 
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4.0  CONCLUSIONS 

Before  restating  the  main  conclusions  of  this  work,  a  few  words  are  in  order 
regarding  the  sensitivity  of  the  results  with  respect  to  variations  of  the  input 
parameters.  The  principal  uncertainties  concern  the  magnitudes  and  the  temper¬ 
ature  dependence  of  the  recombination  coefficients.  These  magnitudes  affect 
TBLIP’  as  seen  ln  Figure  2.  It  is  claimed  here  only  that  several  orders  of  mag¬ 
nitude  reduction  in  the  recombination  coefficient,  which  is  expected  physically 
for  repulsive  centers,  will  produce  the  dramatic  effect  on  Tgrlp.  The  temperature 
dependence  of  the  recombination  coefficients  mostly  affects  the  shape  and  not  the 
magnitude  of  the  calculated  curves.  The  main  conclusions  of  this  work  are  based 
on  physical  arguments  regarding  carrier  dynamics  and  are  largely  insensitive  to 
the  numerical  inputs. 


The  model  employed  in  this  work  indicates  large  improvement  in  TBLIp  for  the 
counterdoped  system,  on  the  order  of  tens  of  degrees.  The  best  responsivities  and 
detectivities  are  obtained  for  exact  compensation  of  the  deep  donor,  although  the 
responsivities  are  not  very  sensitive  to  the  variation  of  the  compensation  ratio 
around  unity.  It  is  recommended  that  the  shallowest  possible  acceptor  be  used  to 
compensate  the  deep  donor  in  order  to  enhance  the  responsivity . 

5 . 0  ACKNOWLEDGEMENTS 

The  author  is  indebted  to  Dr.  Melvin  C.  Ohmer  (AFWAL/MLPO,  Wright-Patterson 
Air  Force  Base,  GH)  and  Mr.  John  A.  Detrio  (University  of  Dayton  Research  Insti¬ 
tute,  Dayton,  OH)  for  many  useful  discussions  and  for  their  encouragement. 


,  t-'i  ■  A  3  J  5  i 


John  A.  Baker 


Universal  Energy  Systems,  Inc. 
4401  Dayton-Xenia  Road 
Dayton,  Ohio  45432 


May  1991 


Final  Report  for  Period  April  1986  -  April  1987 


Approved  for  public  release;  distribution  is  unlimited. 


MATERIALS  DIRECTORATE 

WRIGHT  LABORATORY 

AIR  FORCE  SYSTEMS  COMMAND 

WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433-6533 


/  Unclassified 


SECURITY  CLASSIFICATION  OF  TmIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 


Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


Form  Approved 
OMB  No.  0704-0138 


1b.  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  Public  Release; 
Distribution  is  Unlimited 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

WL-I'R-  91-4048 


6a.  NAME  OF  PERFORMING  ORGANIZATION  I  6b.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 

Universal  Energy  Systems,  IncJ  Of  applicable)  Materials  directorate  (wL/MLPu) 


Universal  Energy  Systems,  IncJ 


6c.  ADDRESS  (City,  State,  and  ZIP  Code) 

4401  Dayton-Xenia  Rd 
Dayton  UH  454311 


8a.  NAME  OF  FUNDING  /  SPONSORING 
ORGANIZATION 

Materials  Directorate 


Sc.  ADDRESS  (City,  State,  and  ZIP  Code) 


Wright  Laboratory 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 

Wri gnt-Patterson  mFB  uH  45433-6533 


8b.  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

WL/MLPU  F33615-86-C-5091 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

Wright  Laboratory 

ELEMENT  NO. 

NO 

NO 

ACCESSION  NC 

Wri ght-Patterson  AFB  OH  454j3-b533 

65502F 

30U5 

50 

79 

11.  TITLE  (Include  Security  Classification) 

Counterdoped  High  temperature,  Silicon  Array  Infrared  Detector 


12.  PERSONAL  AUTHOR(S) 

Baker,  John  A, 


13a.  TYPE  OF  REPORT 

Final  Report 


16.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED 
FROM  4/86  TO. 


14.  DATE  OF  REPORT  (Year,  Month,  Day)  15.  PAGE  COUNT 

Mav  1991  83 


FIELD 

GROUP 

17 

!  05 

COSATI  CODES 


SUB-GROUP 


18  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Counterdoping,  silicon,  electron  irradiation,  infrared 
.  detector,  extrinsic 


19  ABSTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

A  counterdoped  infrared  detector  system  consisting  of  a  silicon  substrate  doped  with  boron 
and  divacancies  has  been  created. 


High  purity  polycrystal  1 i ne  silicon  was  further  purified  by  float  zoning  to  reduce  the 
impurity  levels  still  present;  it  was  then  back  doped  with  boron.  Llectron  irradiation  was 
used  to  create  the  divacancies  which  become  the  active  infrared  centers.  Boron  served  to 
compensate  the  newly  created  di vacancies  thus  producing  only  positively  charged  di vacancies 
these  positively  cnarged  divacancies  are  stable  under  illumination.  The  experimental 
results  thus  verify  the  theoretical  model  postulated  by  Elliott  for  counterdoping. 

A  theoretical  analysis  of  multivalent  statistics  as  applied  to  counterdoped  detectors  has 
developed  an  extended  model  for  counterdoped  detectors.  This  model  has  identified  the 
relevant  parameters  which  control  the  operating  temperature  and  figure  of  merit  for  the 


20  DISTRIBUTION /  AVAILABILITY  OF  ABSTRACT 
ED  UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Gail  J.  Brown 


21  ABSTRACT  SECURITY  CLASSIFICATION 

C3  dtic  users  Unclassified 


22b  TELEPHONE  (Include  Area  Code)  22 c  OFFICE  SYMBOL 

(513)  355-4474  WL/MLPO 


DO  Form  1473,  JUN  86 


Previous  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

Unci assi tied 


block  ly. 

counterdoped  detector.  The  model  may  be  used  to  optimize  doping  concentration  for  the 
highest  background  limited  infrared  photoconductor  (BLIPJ  operating  temperature. 


